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Abstract 

These lectures provide a pedagogical review of the present status of theories explaining 
the observed baryon asymmetry of the Universe. Particular emphasis is given on GUT 
baryogenesis and electroweak baryogenesis. The key issues, the unresolved problems 
and the very recent developments, such as GUT baryogenesis during preheating, are 
explained. Some exercises (and their solution) are also provided. 
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1 Introduction 



We do not know the history of the observable Universe before the epoch of nucleosynthe- 
sis, but it is widely believed that there was an early era of cosmological inflation. The 
attraction of this paradigm is that it can set the initial conditions for the subsequent hot 
big-bang, which otherwise have to be imposed by hand. One of these is that there be no 
unwanted relics (particles or topological defects which survive to the present and contradict 
observation). Another is that the initial density parameter should have the value £1 = 1 to 
very high accuracy, to ensure that its present value has at least roughly this value. There 
is also the requirement that the Universe be homogeneous and isotropic to high accuracy. 
The flatness and the horizon problems of the standard big bang cosmology are - indeed - 
elegantly solved if during the evolution of the early Universe the energy density happens to 
be dominated by the vacuum energy of a scalar field - the inflaton - and comoving scales 
grow quasi-exponentially. 

At the end of inflation the energy density of the Universe is locked up in a combination of 
kinetic energy and potential energy of the inflaton field, with the bulk of the inflaton energy 
density in the zero-momentum mode of the field. Thus, the Universe at the end of inflation 
is in a cold, low-entropy state with few degrees of freedom, very much unlike the present 
hot, high-entropy universe. The process by which the inflaton energy density is converted 
into radiation is known as reheating. What is crucial about these considerations is that, at 
the end of inflation, the Universe does not contain any matter and - even more important 

- the Universe looks perfectly baryon symmetric - there is no dominance of matter over 
antimatter. 

The observed Universe - however - is drastically different. We do not observe any 
bodies of antimatter around us within the solar system and if domains of antimatter exist 
in the Universe, they are separated from us on scales certainly larger than the Virgo cluster 
(~ 10 Mpc). The Universe looks baryon asymmetric to us. Considerations about how the 
light element abundances were formed when the Universe was about 1 MeV hot lead us to 
conclude that the difference between the number density of baryons and that of antibaryons 
is about 1CT 10 if normalized to the entropy density of the Universe. 

Theories that explain how to produce such a tiny number go generically under the name 
of Theories of baryogenesis and they represent perhaps the best example of the perfect 
interplay between particle physics and cosmology Until now, many mechanisms for the 
generation of the baryon asymmetry have been proposed and we have no idea which is 
the correct one. Grand Unified Theories (GUTs) unify the strong and the electroweak 
interactions and predict baryon number violation at the tree level. They are - therefore 

- perfect candidates for a theory of baryogenesis. There, the out-of-equilibrium decay 
of superheavy particles can explain the observed baryon asymmetry, even though there 
remain problems strictly related to the dynamics of reheating after inflation. In the theory 
of electroweak baryogenesis, baryon number violation takes place at the quantum level due 
to the chiral anomaly. Baryogenesis scenarios at the electroweak scale have been the subject 
of intense activity in the last few years. They are certainly attractive because they can be 
tested at the current and future accelerator experiments. 

The bottom line of all this intense research is that, within the standard model of weak 
interactions, it is difficult, if not impossible, to explain how the generation of the baryon 
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asymmetry took place. Therefore, the observation of a baryon asymmetry in the Universe 
is an indication that the description of Nature cannot be limited to the Weinberg-Salam 
theory, something else is called for. 

The goal of these lectures is to provide a pedagogical review of the present state of 
baryogenesis, with particular emphasis on GUT baryogenesis and electroweak baryogenesis. 
The technical details of the numerous models considered in the literature are not elaborated, 
but the key points, the unresolved problems and the very recent developments - such as 
GUT baryogenesis during preheating - are presented. We hope that this approach will help 
the reader to get interested in this fascinating subject. A different focus may be found in 
other accounts of the subject [29, 113, 122, ^q| . Some exercises (and their solution) are also 
provided. 

The review is laid out as follows. Section 2 describes some necessary tools of equilibrium 
thermodynamics. Section 3 contains some considerations about the baryon symmetric Uni- 
verse and explains the three basic conditions necessary to generate the baryon asymmetry. 
The standard out-of-equilibrium scenario is addressed in Section 4, while Section 5 con- 
tains informations about GUT baryogenesis and the thermal history of the Universe, with 
particular attention paid to the recent developments related to the theory of preheating. 
Section 6 is dedicated to the issue of baryon number violation in the standard model and its 
possible implications for GUT baryogenesis and leptogenesis. Finally, section 7 addresses 
the rapidly moving subject of electroweak baryogenesis. 

A note about conventions. We employ units such that % = c = k = 1 and references are 
listed in alphabetic order. 



2 Some necessary notions of equilibrium thermodynamics 

2.1 Expansion rate, number density, and entropy 

Before launching ourselves into the issue of baryon asymmetry production in the early 
Universe, let us just remind the reader a few notions about thermodynamics in an expanding 
Universe that will turn out to be useful in the following. According to general relativity, the 
space-time evolution is determined via the Einstein equation by the matter content of the 
Universe, which differs from epoch to epoch depending on what kind of energy dominates the 
energy density of the Universe at that time. There are three important epochs characterized 
by different relation between the energy density p and the pressure p: 1) vacuum energy 
dominance with p = —p, 2) massless (relativistic) particle dominance with p = p/3 and 3) 
nonrelativistic particle dominance with p <C p- The Einstein equation reads 

Rfw ~ = SttGnT^, (1) 

where R^ v is the Ricci tensor, 7Z is the Ricci scalar, g^ v is the metric, Gn = Mp 2 = 
(1.2 x 10 19 ) -2 GeV -2 is the Newton constant and Tu V is the stress-energy tensor. 

With the homogeneity and isotropy of the three-space, the Einstein equation is much 
simplified with the Robertson- Walker metric 

ds 2 = dt 2 - a 2 (t)x 2 , (2) 
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where a(t) is the cosmic scale factor and the stress-energy tensor is reduced to Tn V = 
~V9iiv + {p + p)UpU v . Here is the velocity vector which in the rest frame of the plasma 
reads = (1,0) and has the property u^u^ = 1. The — component of eq. (||) becomes 
the so-called Friedmann equation 

2 k 8ttG n 
H +^ = —P, (3) 

where k can be chosen to be +1, —1 or for spaces of constant positive, negative or zero 
spatial curvature, respectively, and we have defined the Hubble parameter 

H S *, (4) 
a 

which meaures how fast the Universe is expanding during the different stages of its evolution. 

The fi = component of the conservation of the stress-energy tensor (Tfj? = 0) gives 
the first law of thermodynamics in the familiar form 

d(pa 3 ) = -pd{a 3 ), (5) 

that is, the change in energy in a comiving volume element, d(pa 3 ) is equal to minus the 
pressure times the change in volume, p d(a 3 ). For a simple equation of state p = wp, 
where w is independent of time, the energy density evolves like p oc a~ 3 ^ l+w \ Examples of 
interest include radiation (p oc a~ 4 ), matter (p oc a~ 3 ), vacuum energy (p oc constant). The 
time-behaviour of the scale factor a(t) then is 



a oc e m , H 



2) a t 1 / 2 , 

3) a oc t 2/3 . 



I8itV 
3M|' 



(6) 



The first stage is the inflationary epoch where the constant vacuum energy V gives the 
exponential growth of the scale factor, which is believed to solve the horizon and the flatness 



problems of the standard big-bang theory of cosmology [54] (for a review, see [86]). Of great 
importance is the transient stage from inflation to radiation dominance. This epoch is called 
reheating after inflation and we shall come back to it later in these lectures. 

What is relevant for us is that the early Universe was to a good approximation in thermal 
equilibrium at temperature T [73] and we can define the equilibrium number density rtv 



of a generic interacting species X as 

n x Q = / /eq(p,Mx) d3 P> ( 7 ) 

where gx denotes the number of degrees of freedom of the species X and the phase space 
occupancy /eq is given by the familiar Fermi-Dirac or Bose-Einstein distributions 

/eq(p, nx) = [eM(E x - fix)/?) ± 1] , (8) 
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where Ex = (p 2 + is the energy, px is the chemical potential of the species and 

+1 pertains to the Fermi-Dirac species and —1 to the Bose-Einstein species. 
In the relativistic regime T ^> mx , px formula (0) reduces to 

eq_J (C(3)/7T 2 )g x T 3 (Bose), 
nx \ (3/4)(C(3)/^ 2 ) 5x T 3 (Fermi), W 

where £(3) — 1.2 is the Riemann function of 3. In the non-relativitic limit, T <C mx, the 
number density is the same for Bose and Fermi species and reads 



nT = 9x{^P) 12 e~- + -. (10) 

It is also important to define the number density of particles minus the number density of 
antiparticles 



eq eq 9x 



n x -n T 



(2vr) : 



/ /eq(P,A*x) d 3 p - {fix -fix) 

"gj? [vr 2 + (^) 3 ] (T » mx) 

2g x (m x T/2TTf/ 2 smh(p x /T) exp(-m x /T) (T « m x ). 



(11) 



Notice that, in the relativistic limit T S> mx, this difference scales linearly for T > px- 
This means that detailed balances among particle number asymmetries may be expressed 
in terms of linear equations in the chemical potentials. 

We can similarly define the equilibrium energy density p x of a species X as 



pI q = ^3 / t\r m ( P )<f>,>. (i2) 



which reads in the relativistic limit 

EQ f (vr 2 /30) 5x T 4 (Bose), 
Px \ (7/8)(vr 2 /30)<7xT 4 (Fermi). 1 j 

Since the energy density of a non-relativistic particle species is exponentially smaller than 
that of a relativisitic species, it is a very convenient approximation to include only relativistic 
species with energy density pr in the total energy density p of the Universe at temperature 
T 2 

P^ PR = ^g*T\ (14) 
where g* counts the total number of effectively massless degrees of freedom of the plasma 

9*= E^(^) 4 + ^E^(^) 4 - as) 

i=bos i=ier 

Here Tj denotes the effective temperature of any species % (which might be decoupled from 
the thermal bath at temperature T). In the rest of these lectures we will be always concern 
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with temperatures higher than about 100 GeV. At these temperatures, all the degrees of 
freedom of the standard model are in equilibrium and g* is at least equal to 106.75. 

From this expression we derive that, when the energy density of the Universe was dom- 
inated by a gas of relativistic particles, p oc a -4 oc T 4 and, therefore 73] 



Tex a" 1 . (16) 

Assuming that during the early radiation-dominated epoch (t < 4x 10 10 sec), the scale factor 
scales like t a , where a is a constant, the Hubble parameter scales like oc T 2 oc a~ 2 . This 
means that the scale factor a(t) scales like t 1 / 2 and we recover 2) of Eq. (0). More precisely, 



the expansion rate of the Universe H is [73] 



h=(^op] 1/2 ^i-^gl /2 ^- (17) 



Using the fact that H = (l/2i) and Eq. (|l7|), we can easily relate time and temperature as 

t * 0-301_^- * ™ sec - ( 18 ) 



jV 2 T 2 VMeV 

Another quantity that will turn out to be useful in the following is the entropy density. 
Throughout most of the history of the Universe, local thermal equilibrium is attained and 
the entropy in a comoving volume element s remains constant. Since it is dominated by the 
contribution of relativisitic particles, to a very good approximation 

s = ^r9*sT 3 , (19) 



where 

9*s= E wf^y + sEftff ) • ( 2 °) 



a. I , ,., , 

v T J 8 .^f 36 V T 



For most of the history of the Universe, however, all the particles have the same temperature 
and we can safely replace g*s with g*. Notice that the conservation of entropy implies that 
s oc a -3 and therefore g^sT^a 3 remains a constant as the Universe expands. This means 
that the number of some species X in a comoving volume Nx = a?nx is proportional to 
the number density of that species divided by s, Nx oc nx/ s. 

2.2 Local thermal equilibrium and chemical equilibrium 

So far we have been using the fact that, throughout most of the history of the Universe, 
thermal equilibrium was attained. The characteristic time tx for particles of a species X 
with respect to the process X + A ■ ■ ■ — * C + D + ■ ■ ■ is defined by the rate of change of the 
number of particles per unit volume nx due to this process 

±—±(*2L) . (21) 

tx nx \ at ) x+a- — >c+d+- 
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In the early Universe, if tx is smaller than the characteristic time of the expansion H -1 , 
then there is enough time for the process to occur and the particles X's are said to be 
thermally coupled to the cosmic fluid. By contrast, if tx 3> H -1 , for every process in which 
the particles X's are involved, then they are not in thermal equilibrium and they are said 
to be decoupled. 

In order to analyze the evolution of the particle populations which constitute the cosmic 
fluid, it is necessary to compare H~ l with tx at different temperatures. This is done through 
the Boltzmann equation ||], which, in an expanding Universe, reads 

' d 'a 3 n x )= [ n x C[fx], (22) 



a 3 dt 

where C is the collision operator. Eq. ( p2|) may be rewritten as 

^ + 3Hn x = £ J nftf m — {l±fx)0-±fj)~'W(e + m + — ^X+j + .~) 

- fxf j ---(l±h)(l±f m )---W(X + j + ---^e + m + ---), (23) 

where ir = nx^j • ■ ■ ^e^m • • •, 7Tj = (2ir)~ 3 gi(d 3 p/2Ei) is the volume element in the phase 
space, W is the matrix element of the given process and (+) applies to bosons and (— ) to 
fermions.. The second term in the left-hand side of Eq. (|23|) accounts for the nx diluition 
due to the cosmic expansion and the right-hand side accounts for the nx variations due to 
any elemenatry process X + j + • • • — > £ + m + ■ ■ ■ in which the X particles are involved. 
As it stands, Eq. ( |2~3| ) is rather formidable and complicated, but some approximations can 
be made to transform it in a simpler form. 

Let us consider, for example, a process like X + / — > X' + /, where the number of X 
particles does change in the scatterings and let us also suppose that the / particles are light 
(T > mj) and that the corresponding population is in thermal equilibrium. In the case in 
which the X distribution function is described by a Maxwell-Boltzmann distribution, i.e. 
the X particles are in equilibrium at temperatures smaller than mx, it is easy to see that 



the right-hand side of Eq. (23) may be expressed in the form 

r.h.s. of Eq. (23) = -\n x - n x Q ] S (24) 

where 



S = J n f E f ff Q a(X + f^X + f) 



n 



EQ 

/ 



(a(X + f^X' + f)v). (25) 



The notation (av) stands for the thermal average cross section times the relative velocity 
v. The inverse time scale t x 1 = Tx associated to the elastic process is therefore 

Tx = r x 1 ^n^ Q (a(X + f^X f + f)v). (26) 

From these very simple considerations, we may conclude that the X degrees of freedom are 
in thermal equilibrium if 

T x ^ n^ Q (a(X + / X + f)v) > H (thermal equilibrium is attained). (27) 
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Departure from thermal equilibrium is expected whenever a rate crucial for mantaining 
thermal equilibrium becomes smaller than the expansion rate, Tx < H. 

Another useful concept is that of chemical equilibrium. In general, a species X is in 
chemical equilibrium if the inelastic scatterings which change the number of X particles 
in the plasma, X + j — > t + m, have a rate rj ne i larger than the expansion rate of the 
Universe. In such a case, one is allowed to write down a relation between the different 
chemical potentials /x's 

fj>x + y-j = fJ>e + (28) 

of the particles involved in the process. With these simple notions in mind we may start 
our voyage towards the country of baryogenesis. 

3 The graveyard for a baryon symmetric Universe 

The CPT theorem assures that any particle species X there exists the antiparticle X with 
exactly the same mass, mx = rrl x~- l and decay width, Tx = T-^, and eventually opposite 
charges associated to these particles, Qx = —Q~x- This striking symmetry would naturally 
lead us to conclude that the Universe contains particles and antiparticles in equal number 
densities, nx = n^- The observed Universe, however, is drastically different. We do not 
observe any bodies of antimatter within the solar system and only antiprotons p in the 
cosmic rays, which are believed to be of extra solar origin. Antiprotons are likely to be 
produced as secondaries in collisions pp — > 3p + p at a rate similar to the observed one 

^ ~ 3 x 10" 4 . (29) 
n p 

The experimental limit on riin e /n^ e is similarly of the order of 10~ 5 . We cannot exclude, of 
course, that the dominance of matter over antimatter is only local and is only realized up to 
a certain length scale £b, beyond which the picture is reversed and islands of antimatter are 
found. However, the size of our matter domain must be quite large, roughly speaking > 
10 Mpc [|118| , |119f| (for more restrictive bounds see |^6|). Indeed, for smaller scales one would 
expect a significant amount of energetic 7-rays coming from the reaction of annihilation of 
pp into 7r-mesons followed by the subsequent decay 7r° — ► 27, which would take in the 
boundary area separating the matter and antimatter islands. Another signature for the 
presence of domains of antimatter would be the distortion of the spectrum of the cosmic 
microwave background radiation. In such a case, the permitted value of £b might be smaller 
if voids separate matter and antimatter domains. These voids might be created because of 
an excessive pressure produced by the annihilations at earlies stages of the evolution of the 
Universe or because of low density matter and antimatter in the boundary regions, provided 
that the baryon asymmetry changes sign locally so that in the boundaries it is zero or very 
small. 

All these considerations lead us to conclude that, if domains of matter and antimatter 
exist in the Universe, they are separated on scales certainly larger than the radius of our 
own galaxy (~ 3 Kpc) and most probably on scales larger than the Virgo cluster (~ 10 
Mpc). A much more severe bound on £b (~ 300 Mpc) is potentially reachable by the Alpha 
Magnetic Spectrometer (AMS) @, a detector for extraterrestrial study of antimatter, matter 
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and missing matter which, after a precursor flight on STS91 in May 1998, will be installed 
on the International Space Station where it should operate for three years. 



3.1 Some considerations on nucleosynthesis and the baryon number 

The baryon number density does not keep constant during the evolution of the Universe 
because it scales like a -3 , where a is the cosmological scale factor [73|. It is therefore 
convinient to define the baryon asymmetry of the Universe in terms of the quantity 

V = — , (30) 

n 7 

where ns = — is the difference between the number of baryons and antibaryons per 
unit volume and n 7 = 2^P-T 3 is the photon number density at a certain temperature T. 
The parameter r\ is essential for determining the present light element abundances produced 
at the nucleosynthesis epoch. The parameter r/ may have not changed since nucleosynthesis. 
At these energy scales (~ 1 MeV) the baryon number is conserved if there are no processes 
which would have produced entropy to change the photon number. 
Let us now estimate rj. The baryon number density is 

UB = = Pc, (31) 

m B m B 

where ps is the baryonic energy density and Qb = Pb/ Pc- Using the critical density 

p c = 1.88 x 1(T 29 h 2 gr cm" 3 , (32) 

where 0.5 <h< 0.9 parametrizes the present value of the Hubble parameter Hq, h = H/100 
Km Mpc -1 sec , we obtain 

n B = l.lxlO _5 /i 2 n B cm" 3 . (33) 

On the other hand, the present temperature of the background radiation is To = 2.735 °K 
giving rise to 

** m (mbg)' aa ~*- <34> 



Putting (33) and ( p3|) together, we obtain 

Tr x ~ 3 



rj = 2.65 x 10- 8 n B h 2 ( 2 7 ^ QR ) . (35) 



The range of rj consistent with the deuterium and 3 He primordial abundances is [73] 

4(3) x 1(T 10 < 77 < 7(10) x 1(T 10 , (36) 

where the most conservative bounds are in parenthesis. Conversely we may write the range 
for Qb h 2 to be 

0.015(0.011) < n B h 2 < 0.026(0.038). (37) 
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Sometimes it is useful to describe the baryon asymmetry in terms of B = tlb/s, where s is 
the entropy density of the Universe at a certain temperature T. The range ( |36|) translates 
into 

5.7(4.3) x KT 11 <B< 9.9(14) x HT 11 . (38) 

Now, the fundamental question is: are we able to explain the tiny value of r\ within the 
standard cosmological model? 

Suppose that initially we start with r] = 0. We can compute the final number density 
of nucleons b that are left over after annihilations have frozen out. At temperatures T < 1 
GeV the equilibrium abundance of nucleons and antinucleons is 

-^-^Uf) e-- (39) 

When the Universe cools off, the number of nucleons and antinucleons decreases as long 
as the annihilation rate Tann 

~ nb{aAv) is larger than the expansion rate of the Universe 
H ~ 1.66 gl^ 2 'j^-- The thermally averaged annihilation cross section (<jav) is of the order 
of ml . At T ~ 20 MeV, 

Tann — H and annihilations freeze out, nucleons and antinucleons 
being so rare that they cannot annihilate any longer. Therefore, from (|39|) we obtain 

~ m ~ io-i8 ; (40) 

which is much smaller than the value required by nucleosynthesis. In order to avoid the 
annihilation catastrophe, we may suppose that hypothetical new interactions separated 
matter from antimatter before T ~ 38 MeV, when r] ~ 10~ 10 . At that time, t ~ 10 -3 sec, 
however, the causal region (horizon) was small and contained only ~ 1O _7 M0. Hence we 
cannot explain the asymmetry over the galaxy scales. This argument is not valid, however, 
in cosmological models invoking inflation. Indeed, in these models the region of the Universe 
which is causally connected today was connected even at times ~ 10 -3 sec. These scenarios 
pose other serious cosmological drawbacks, though. If the processes responsible for the 
separation of matter from antimatter took place before inflation, then the baryon number 
was diluted by an enormous factor ~ exp(200), because of the entropy production due to 
inflation. On the other side, if the separation took place after inflation, then it is not clear 
how to eliminate the boundaries separating matter from antimatter islands. 

Another possibility may be represented by explaining the tiny value of rj via statistical 
fluctuations in the baryon and antibaryon distributions. Our own galaxy containes at the 
present epoch approximately 10 79 photons. The comoving volume V that encompasses our 
galaxy today contains about 10 69 baryons, but when the temperature was T > 1 GeV, it 
contained about 10 79 baryons and antibaryons. Frome pure statistical fluctuations one may 
expect an asymmetry (n& — rer)/n& ~ (n^U) -1 / 2 ~ 10 -39 ' 5 , which is again far too small to 
explain the observed baryon asymmetry. 

In conclusion, in the standard cosmological model there is no explanation for the small- 



ness of the ratio (36), if we start from rj = 0. An initial asymmetry may be imposed by 
hand as an initial condition, but this would violate any naturalness principle and would be 
extremely boring! 
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3.2 The three basic conditions for baryogenesis 



As we have already learned, the Universe was initially baryon symmetric (rib = rir) although 
the matter-antimatter asymmetry appears to be large today (rib 3> ni). In the standard 
cosmological model there is no explanation for such a small value of the baryon asymmetry 
consistent with nucleosynthesis and it has to be imposed by hand as an initial condition. 
This option is far from being appealing. However, it has been suggested by Sakharov long 
ago [114] that a tiny baryon asymmetry B may have been produced in the early Universe. 
Three are the necessary conditions for this to happen. 



Exercise 1 



Show that the baryon asymmetry is zero if there is no baryon number violation. 



3.2.1 Baryon number violation 

This condition is somehow obvious since we want to start from a baryon symmetric Universe 
(B = 0) and to evolve it to a Universe where B / 0. Baryon number violation interactions 
are therefore mandatory. They might also mediate proton decay; in such a case phenomeno- 
logical constraints are provided by the lower bound on the proton lifetime r p > 5 x 10 32 
years.. 



3.2.2 C and CP violation 

C (charge conjugation symmetry) and CP (the product of charge conjugation and parity) 
are not exact symmetries. Indeed, were C an exact symmetry, the probability of the process 
i —* f would be equal to the one of the process i — > f. Since the baryon number of / is 
equal in absolute value and opposite in sign to that of /, the net baryon number B would 
vanish. C is maximally violated by the weak interactions. 

Furthermore, because of the CPT theorem, CP invariance is equivalent to time-invariance 
(time reversal). The latter assures that the rate of the process 

i(n, Pi, -> f(rj, Pj,Sj) (41) 

and that of its time-reversed process 

f(rj, ~Pj, -Sj) -> i(r i3 -pi, -s^ (42) 

are equal. Thus, even though it is possible to create a baryon asymmetry in a certain 
region of the phase space, integrating over all momenta p and summing over all spins s 
would produce a vanishing baryon asymmetry. CP violation has been observed in the kaon 
system. However, a fundamental understanding of CP violation is still lacking. Hopefully, 
studies of baryogenesis may shed some light on it. 



3.2.3 Departure from thermal equilibrium 

If all the particles in the Universe remained in thermal equilibrium, then no preferred 
direction for time may be defined and the CPT invariance would prevent the appearance 
of any baryon excess, making the presence of CP violating interactions irrelevant. 
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Let us suppose that a certain species X with mass mx is in thermal equilibrium at 
temperatures T <C mx- Its number density will be given by 

o /o m X V 1 X 

n x ~ gx{m x Tfl 2 e~— +— , (43) 

where fj,x is the associated chemical potential. 

As we have mentioned in the previous Section, a species X is in chemical equilibrium if 
the inelastic scatterings which change the number of X particles in the plasma, X + A — > 
B + C, have a rate ri ne i larger than the expansion rate of the Universe. In such a case, one 
can write down a relation among the different chemical potentials of the particles involved 
in the process 

Vx + = fJ>B + MC- (44) 
In this way the number density in thermal equilibrium of the antiparticle X (mx = m x ) is 

o In m X ftX 

n x — 9x(mxT) ' e t ~, (45) 
where we have made use of the fact that pL x = — fix because of the process 

XX -» 77, (46) 
and /i-y = 0. If the X particle carries baryon number, then B will get a contribution from 

B oc n x - n x = 2g x (m x T) 3 / 2 e-^smh (^f) . (47) 

The crucial point is now that, if X and X undergo i?-violating reactions, as required by 
the first Sakharov condition, 

XX -> X~X, (48) 

then fi x = and the relative contribution of the X particles to the net baryon number 
vanishes. Only a departure from thermal equilibrium can allow for a finite baryon excess. 



4 The standard out-of-equilibrium decay scenario 

Out of the three Sakharov conditions that we discussed in the previous section, the baryon 
number violation and C and CP violation may be investigated thoroughlly only within 
a given particle physics model, while the third condition - the departure from thermal 
equilibrium -may be discussed in a more general way. Very roughly speaking, the various 
models of baryogenesis that have been proposed so far fall into two categories: 

- models where the out-of-equilibrium condition is attained thanks to the expansion of 
the Universe and the presence of heavy decaying particles; 

- models where the departure from thermal equilibrium is attained during the phase 
transitions which lead to the breaking of some global and/or gauge symmetry. 

In this lecture we will analyse the first category -the standard out-of-equilibrium decay 
scenario {73|, 
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4.1 The conditions for the out-of-equilibrium decay scenario 

It is obvious that in a static Universe any particle, even very weakly interacting, will attain 
sooner or later thermodynamical equilibrium with the surroinding plasma. The expansion 
of the Universe, however, introduces a finite time-scale, tjj ~ H . Let suppose that X 
is a baryon number violating superheavy boson field (vector or scalar) which is coupled to 

1/2 

lighter fermionic degrees of freedom with a strength a x (either a gauge coupling a gaU g e or 
a Yukawa coupling ay). 

In the case in which the couplings are renormalizable, the decay rate Tx of the super- 
heavy boson may be easily estimated to be 

T x ~a x M x , (49) 

where Mx is the mass of the particle X. In the opposite case in which the boson is a gauge 
singlet scalar field and it only couples to light matter through gravitational interactions - 



this is the case of singlets in the hidden sector of supergravity models []103[| - the decay rate 
is from dimensional arguments 

rv ~ £|. (50) 

At very large temperatures T 3> Mx, it is assumed that all the particles species are in 
thermal equilibrium, i.e. nx — n-^ ~ n 7 (up to statistical factors) and that B = 0. At 
T < Mx the equilibrium abundance of X and X relative to photons is given by 



e-~ , (51) 



where we have neglected the chemical potential nx ■ 

For the X and X particles to mantain their equilibrium abundances, they must be able 
to diminish their number rapidly with respect to the Hubble rate H(T). The conditions 
necessary for doing so are easily quantified. The superheavy X and X particles may attain 
equilibrium through decays with rate Tx, inverse decays with rate T 1 ® 

r iD_ r / 1 T>M X , ( . 

x ~ X \ (M X /T)3/ 2 exp(-M x /T) T < M x , ^ 

and annihilation processes with rate T^p 11 oc nx- The latter, however are "self-quenching" 
and therefore less important than the decay and inverse decay processes. They will be 
ignored from now on. Of crucial interest are the i?-nonconserving scattering processes 
2 «-> 2 mediated by the X and X particles with rate T x 

^"^"^ (Mi+W <53> 



where a ~ <? 2 /47r denotes the coupling strength of the X boson. At high temperatures, the 

4- 



2 <-> 2 scatterings cross section is a — a 2 /T 2 , while at low temperatures a — a 2 T 2 jM\ 
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For baryogenesis, the most important rate is the decay rate, as decays (and inverse 
decays) are the mechanism that regulates the number of X and X particles in the plasma. 
It is therefore useful to define the following quantity 

(54) 

T=M X 

which measures the effectiveness of decays at the crucial epoch (T ~ Mj) when the X and 
X particles must decrease in number if they are to stay in equilibrium. Note also that for 
T < Mx, K determines the effectiveness of inverse decays and 2^2 scatterings as well: 
T^/H ~ (Mx/Tfl 2 exp(-M x /T) K and T%/H ~ a{T/M x f K. 
Now, if K ^> 1, and therefore 

Tx » H\ T=Mx , (55) 

then the X and X particles will adjust their abundances by decaying to their equilibrium 
abundances and no baryogenesis can be induced by their decays -this is simply because 
out-of-equilibrium conditions are not attained. Given the expression (|17]) for the expansion 
rate of the Universe, the condition (|55|) is equivalent to 

M x < 5*~ 1/2 a x M P (56) 
for strongly coupled scalar bosons, and to 

M x » gl' 2 M P , (57) 

for gravitationally coupled X particles. Obviously, this last condition is never satisfied for 
M x <M P . 

However, if the decay rate is such that Kcl, and therefore 

Tx < H\ T=Mx , (58) 

then the X and X particles cannot decay on the expansion time-scale tjj and so they remain 
as abundant as photons for T < Mx- In other words, at some temperature T > Mx, the 
superheavy bosons X and X are so weakly interacting that they cannot catch up with the 
expansion of the Universe and they decouple from the thermal bath when still relativistic, 
nx — nj ~ ra 7 at the time of decoupling. Therefore, at temperature T ~ Mx, they will 
populate the Universe with an abundance which is much larger than the equilibrium one. 
This overbundance with respect to the equilibrium abundance is precisely the departure from 
thermal equilibrium needed to produce a final nonvanishing baryon asymmetry. Condition 
(|58|) is equivalent to 

M x >g* 1/2 a x Mp (59) 
for strongly coupled scalar bosons, and to 

M x <gl /2 Mp, (60) 

for gravitationally coupled X particles. It is clear that this last condition is always satisfied, 
whereas the condition ( |59| ) is based on the smallness of the quantity g* ax- In particular, 



K = 



H 
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if the X particle is a gauge boson, ax ~ Ogauge can span the range (2.5xl0 _2 -10 _1 ), while 
is about 10 2 . In this way we obtain from ( |59|) that the condition of out-of-equilibrium 
can be satisfied for 

Mx > (10" 4 - 10" 3 ) M P ~ (10 15 - 10 16 ) GeV. (61) 

If X is a scalar boson, its coupling ay to fermions / with mass is proportional to the 
squared mass of the fermions 

( m f \ 2 
\mw / 

where my/ is the VF-boson mass and ay is typically in the range (10 -2 — 10 -7 ), from where 

M x > (10~ 8 - 1(T 3 ) M P ~ (10 10 - 10 16 ) GeV. (63) 

Obviously, condition dHTJ) is more easily satisfied than condition (|6l| ) and we conclude that 
baryogenesis is more easily produced through the decay of superheavy scalar bosons. On the 
other hand, as we have seen above, the condition (^) tells us that the out-of-equilibrium 
condition is automatically satisfied for gravitationally interacting particles. 



4.2 The production of the baryon asymmetry 

Let us now follow the subsequent evolution of the X and X particles. When the Universe 
becomes as old as the lifetime of these particles, t ~ H^ 1 ~ , they start decaying. This 
takes place at a temperature Tjj defined by the condition 

Tx * H\ t=Td , (64) 

i.e. at 

T D ~ g~ 1/4 o% 2 (MxMp) 1 / 2 < M x , (65) 

where the last inequality comes from (|59|) and is valid for particles with unsuppressed 
couplings. For particles with only gravitational interactions 

1 /2 

T D ~g^ A M x <M X , (66) 

the last inequality coming from (p0[). At T ~ To, X and X particles start to decay and 
their number decrease. If their decay violate the baryon number, they will generate a net 
baryon number per decay. 

Suppose now that the X particle may decay into two channels, let us denote them by a 
and b, with different baryon numbers B a and B^, respectively. Correspondingly, the decay 
channels of X, a and b, have baryon numbers — B a and —B^, respectively. Let r(r) be the 
branching ratio of the X(X) in channel a(a) and 1 — r(r) the branching ratio of X{X) in 
channel b(b), 

F(X -» a) 
r = —fx—> 
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r(x -» 


a) 






r(x -» 








r(x 







1-r = -i— (67) 



where we have been using the fact that the total decay rates of X and X are equal because 
of the CPT theorem plus unitarity. 

The average net baryon number produced in the X decays is 

rB a + (1 - r)B b , (68) 

and that produced by X decays is 

-rB a -{l-r)B h . (69) 

Finally, the mean net baryon number produced in X and X decays is 

AB = (r — r)B a + [(1 - r) - (1 - r)] B b = (r - r){B a - B b ). (70) 

Equation ( |70|) may be easily generalized to the case in which X(X) may decay into a set of 
final states f n (f n ) with baryon number B n (—B n ) 

AB = -^Y, B n - fn) ~ T(X - /„)] . (71) 

x n 

At the decay temperature, Tjj < Mx, because K <C 1 both inverse decays and 2 «-> 2 baryon 
violating scatterings are impotent and can be safely ignored and thus the net baryon number 
produced per decay AB is not destroyed by the net baryon number —AB produced by the 
inverse decays and by the baryon number violating scatterings. 

At T ~ Tjj, nx — njF — n 7 and therefore the net baryon number density produced by 
the out-of-equilibrium decay is 

n B = ABn x , (72) 

from where we can see that AB coincides with the parameter r\ defined in ([}(]) if nx — n^. 

The three Sakharov ingredients for producing a net baryon asymmetry can be easily 
traced back here: 

- If B is not violated, then B n = and AB = 0. 

-If C and CP are not violated, then T(X — > f n ) = T(X — > f n ), and also AB = 0. 

- In thermal equilibrium, the inverse processes are not suppressed and the net baryon 
number produced by decays will be erased by the inverse decays. 

Since each decay produces a mean net baryon number density n# = ABnx — ABn^ 
and since the entropy density is s ~ g*n^, the net baryon number produced is 

n B ABn 7 AB 
B = — ~ '- ~ . (73) 
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Taking ~ 10 2 , we see that only tiny C and CP violations are required to generate 
AB ~ 1(T 8 , and thus B ~ 1(T 10 . 

To obtain (73) we have assumed that the entropy realese in X decays is negligible. 
However, sometimes, this is not a good approximation (especially if the X particles decay 
very late, at Tjj <C Mx, which is the case of gravitationally interacting particles). In that 
case, assuming that the energy density of the Universe at Td is dominated by X particles 

px^Mxnx, (74) 

and that it is converted entirely into radiation at the reheating temperature Trh 

IT 2 

P=2q9*T rh , (75) 

we obtain 

n * (76> 

We can therefore write the baryon number as 

B ~ ^-^-AB. (77) 
4 M x 

We can relate Trjj with the decay rate Tx using the decay condition 



and so we can write 

/ -1/2 \ 1 / 2 

For the case of strongly decaying particles (through renormalizable interactions) we obtain 
while for the case of weakly decaying particles (through gravitational interactions) we obtain 

(—1/2 \ V 2 

AB. (8!) 

In the other extreme regime K S> 1, one expects the abundance of X and X bosons to 
track the equilibrium values as Tx 3> H for T ~ Mx ■ If the equilibrium is tracked precisely 
enough, there will be no departure from thermal equilibrium and no baryon number may 
evolve. The intermediate regime, K ~ 1, is more interesting and to address it one has 
to invoke numerical analysis involving Boltzmann equations for the evolution of B. This 
has been done in refs. [71, 57, 72]. The numerical analysis essentially confirms the 



qualitative picture we have described so far and its discussion is beyond the scope of these 
lectures. 
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'--<■' 

f 2 

(c) (d) 
Figure 1: Couplings of X and Y to fermions /j. 

4.2.1 An explicit example 

Let us consider first two massive boson fields X and Y coupled to four fermions /i, / 2 , 
fs and through the vertices of Fig. 1 and describing the decays X — > /1/21/3/4 and 
Y -> 73/1,74/2- We will refer to_these_vertices as (f 2 \X\fi), (/4J-XI/3), (/i|y|/ 3 ) and 
(f2\Y\U), and their CP conjugate X -> / 2 /i,/ 4 /3 and y -> f l fs,f 2 U by their complex 
conjugate. In the Born approximation AP» = because from ( ffT| ) one finds 

r(X -> 7l/2)Born = 1$ \ {h \X\ fl ) | 2 = T(X -» 7 2 /l)Born, (82) 

where /j^ 2 accounts for the kinematic structures of the processes X —* fxf% and X — > / 2 /i 
and the same may be found for the other processes contributing to AB. This shows that, 
to obtain a non-zero result for AB, one must include (at least) corrections arising from 
the interference of Born amplitudes of Fig. 1 with the one-loop amplitude of Fig. 2. 
For example, the interference of the diagrams in Fig. 1(a) and Fig. 2(a) (in the square 
amplitude) is shown in Fig. 3(a), where the thick dashed line is the unitarity cut (equivalent 
to say that each cut line represents on-shell mass particles) . The amplitude of the diagram 
in Fig. 3(a) is given by /j s 2 y 4 ^i234, where the kinematic factor ij^y 4 accounts for the 
integration over the final state phase space of fi and / 1 and over momenta of the internal 
states fi and / 3 , and 

^1234 = (/l|y|/3)*(/4|X|/ 3 )(/ 2 |y|/ 4 )(/2|X|/l)*. (83) 
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Figure 2: One-loop corrections to the Born amplitude of Fig. 1. 






Figure 3: Intereference between the diagrams of Fig. 1 and Fig. 2 for the square amplitudes 
of X decay. 
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The complex conjugate diagram of Fig. 3(b) has the complex conjugate amplitude. 
Therefore, the contribution from the diagrams in Figs. 3(a) and 3(b) to the decay X — > f 1 fa 
is 



T(X ► f 1 fa ) interference 



4y 4 ^i234 + h.c. 



(84) 



To obtain the CP conjugate amplitude X — > / 2 /i all couplings must be complex conju- 
gated, although the kinematic factors Ixy are unaffected by CP conjugation. Therefore 
the interference contribution to the X — > f 2 fa decay rate is given by 



T(X - / 2 /x) 

interference 

and the relevant quantity for baryogenesis is given by 

r(x - 7i/ 2 ) - r(X - J 2 /i) = -4 im 



rl234 * , h 
1 XY "1234 + nx - 



1234 
XY 



Im [ft 



1234J 



(85) 



(86) 



The diagrams of the decays X — > / 3 /4 and X — > / 3 /4 differ from the one in Figs. 3(a) and 
3(b) only in that the unitarity cut is taken through fa and fa instead of fa and fa. One 
easily obtains 



F(X - / 3 / 4 ) - r(X - / 4 / 3 ) = -4 Im 



t3412 

'xy 



Im [ft 



1234J 



(87) 



The kinematic factors Ixy for loop diagrams may have an imaginary part whenever any 
internal lines may propagate on their mass shells in the intermediate states, picking the 
pole of the propagator 



PP 



P 



2 — m 2 + ie 



+ m5(p 2 — m 2 ), 



where PP stands for the principal part. This happens if Mj > 7714 +777-2 and > 7773+777,4. 
This means that with light fermions, the imaginary part of Ixy will be always nonzero. The 
kinematic factors Im [ij^y ] an d I m [^xy] are therefore obtained from diagrams involving 
two unitarity cuts: one through the lines fa and fa and the other through the lines fa and 
fa. The resulting quantities are invariant under the interchanges fa <-> fa and fa <-> fa and 
consequently 



Im 



rl234 
1 XY 



Im 



3412 
XY 



Im [J 



xyj 



(89) 



Defining Bi the baryon number of the fermion fa, the net baryon number produced in the 
X decays is therefore 



(AB); 



Im [I XY ] Im [ft 1234 ] [j5 4 - S 3 - (5 2 - Bi)] . 



(90) 



To compute the baryon asymmetry (AB)y one may observe that the set of vertices in Fig. 
1 is invariant under the transformations X <-> Y" and fa ^ fa. These rules yield 



(AB)y = — Im [J yx ] Im [fti" 234 ] [B 4 - B 3 - (#2 - B 1 )} 

1 y 



(91) 
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and the total baryon number is therefore given by 

(AB) = (AB) X + (AB)y = 4| Im i [^ y] - Im ^ yx] | Im [fi 1234 ] [B 4 -B 3 - (B 2 - B 1 )} . 

(92) 

We can notice a few things: 

- If the X and Y couplings were B conserving, the two possible final states in X and Y 
decays would have the same baryon number, i.e. B4 — B3 = B 2 — B\ and therefore AB = 0. 
Therefore the baryon number must be violated not only in X decays but also in the decays 
of the particle exchanged in the loop. 

- Some coupling constants in the Lagrangian must be complex to have Im [01234]- 

- Even if (AB)x and (AB)y are both nonvanishing, the sum can be vanish if the first 
bracket in (^2[) cancels out. This happens if the X and Y particles have the same mass and 
Tx =Ty. 



4.3 Baryon number violation in Grand Unified Theories 

The Grand Unified Theories (for a review, see [^]) try to describe the fundamental in- 
teractions by means of a unique gauge group G which contains the Standard Model (SM) 
gauge group SU (3)c ® SU{2) l®U (l)y . The fundamental idea of GUTs is that at energies 
higher than a certain energy threshold Mgut the group symmetry is G and that, at lower 
energies, the symmetry is broken down to the SM gauge symmetry, possibly through a chain 
of symmetry breakings 

G M ^ G{1) A 4 G{2) SU(3) C ® SU(2) L ® U(l)y, (93) 

corresponding to 

M GU T >M l >M 2 > ■ ■ ■ Mw, G D D G (2) D • • ■ SU(3) C ® SU(2) L ® U(l)y. (94) 

What is the main motivation for invoking GUTs? Gauge couplings (couplings to gauge 
fields) are charactherized by a dimensionless constant g, or equivalently by a = g 2 /Air. (For 
electromagnetism, g is the electron charge and a evaluated at low energy is the fine structure 
constant a em = 1/137.) Gauge couplings are not supposed to be extremely small, and one 
should take g ~ 1 for crude order of magnitude estimates (making a one or two orders of 
magnitude below 1). Assuming small couplings, the perturbative effects usually dominate, 
and we focus on them for the moment. With perturbative quantum effects included, the 
effective masses and couplings depend on the relevant energy scale Q. The dependence on Q 
(called 'running') can be calculated through the renormalization group equations (RGE's), 
and is logarithmic. In the context of collider physics, Q can be taken to be the collision 
energy, if there are no bigger relevant scales (particle masses). For the Standard Model 
there are three gauge couplings, ai where i = 3,2,1, corresponding for respectively to 
the strong interaction (colour SU(3)c) the non-abelian electroweak interaction {SU{2)i) 
and electroweak hypercharge (U(l)y). (The electromagnetic gauge coupling is given by 
a~ l = a^ 1 + a% .) In the one-loop approximation, ignoring the Higgs field, their running 
is given (at one loop) by 

don b { 

:«i • (95) 



dln(Q 2 ) 4vr 
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The coefficients 6j depend on the number of particles with mass <C Q- Including all particles 
in the minimal supersymmetric standard model gives b\ = 11, 62 = 1 an d 63 = —3. 

Using the values of a, measured by collider experiments at a scale Q ~ 100 MeV, one 
finds that all three couplings become equal at a scale |2|, |3^, |82|Q Q = Mgut, where 

M GU T ^ 2 x 10 16 GeV. (96) 

The unified value is 

a G UT ^ 1/25. (97) 

One explanation of this remarkable experimental result may be that there is a GUT, in- 
volving a higher symmetry with a single gauge coupling, which is unbroken above the scale 
Mgut- Another might be that field theory becomes invalid above the unification scale, to 
be replaced by something like weakly coupled string theory or M-theory |124| which is the 



source of unification. At the time of writing there is no consensus about which explanation 
is correct, but in this section we will focus on Grand Unified Theories and their relevance 
for baryogenesis. 

It is a general property of GUTs that the same representation may contain both quarks 
and leptons and therefore there exist gauge bosons which mediate gauge interactions among 
fermions having different baryon number. This is not enough -though- to conclude that 
in GUTs the baryon number is violated, because it might be possible to assign a baryonic 
charge to the gauge bosons in such a way that each vertex boson-fermion-fermion the baryon 
number is conserved. Let us discuss this crucial point in more detail. 

The fundamental fermions of the SM are 

t L = (1,2,-1/2), 

Ql = (3,2,1/6), 

el = (1,1,1), 
u£ = (3,1,-2/3), 

d c L = (3,1,1/3), (98) 

where in parenthesis we have written then SU(3)c 8> SU(2)l <8> U{\)y quantum numbers 
and all the spinors are left-handed. Given two spinors ipL and xl, it is possible to define a 
renormalizable coupling to a gauge boson by 

^[^XL^ + h.c, (99) 

where = (1, a) and a are the Pauli matrices. At this point one may try to write down 



all the couplings of the form (99) starting from the spinors of the SM and identify all the 
possible gauge bosons which may be present in a GUT having the same spinors of the SM. 
Of course, the same gauge boson may be coupled to more than one pair of spinors. If all 
the spinor pairs have the same baryon number B, then it suffices to assign a baryon number 



To be precise, |ai = 012 = Q3 = Qgut, the factor 5/3 arising because the historical definition of 
ai is not very sensible. In passing we note that the unification fails by many standard deviations in the 
absence of supersymmetry, which may be construed as evidence for supersymmetry and anyhow highlights 
the remarkable accuracy of the experiments leading to this result. 
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— B to the gauge boson and obtain a baryon number conserving theory. If there exist gauge 
bosons which couple to spinor pairs having different baryon number, one may write down 
baryon number violating interactions. These bosons are given in Table 1, where we have 
indicated, for every gauge boson, all the possible interactions and the corresponding baryon 
numbers B and baryon minus lepton numbers B — L 



Table 1 



Gauge boson 


spinors 


B 


B-L 


^ = (3,2,-5/6) 




2/3 


2/3 




QVl 


-1/3 


2/3 






-1/3 


2/3 


V 2 = (3,2,1/6) 




-1/3 
2/3 


2/3 
2/3 



Of course, every gauge boson listed in Table 1 has the corresponding antiboson. One 
can repeat the same procedure to identify the scalar bosons 5 which may mediate baryon 
number violation interactions via fermions. The generic coupling reads 

i X T L<J 2 ^ L S + h.c. (100) 

If we consider all the spinor pairs Xl^Pl, even belonging to different families, we get the 
following possibilities 



Table 2 



Scalar boson 


spinors 


B 


B-L 


5 X = (3, 1,-1/3) 




-1/3 


2/3 




QlQl 


2/3 


2/3 




e L u L 


-1/3 


2/3 




44 


2/3 


2/3 


5^ = (3, 1,-4/3) 


e c L d c L 


-1/3 


2/3 




cf cf 
U L U L 


2/3 


2/3 


5 3 = (3,3,-1/3) 




-1/3 


2/3 




QlQl 


2/3 


2/3 



Out of all possible scalar and gauge bosons which may couple to the fermions of the 
SM, only the five that we have listed may give rise to interactions which violate the baryon 
number. A crucial point for what we will be discussing in the following is that each of 
these bosons have the same combination B — L, which means that this combination may 
be not violated in any vertex boson-fermion-fermion. This is quite a striking result and 
originates only from having required the invariance under the SM gauge group and that the 
only fermions of the theory are those of the SM. 

The extension of the fermionic content of the theory may allow the presence of more 
heavy bosons which will possibly violate B and even B — L. In the Grand Unified Theories 
based on 50(10) - for instance - there is another fermion which is a singlet under the SM 
gauge group and is identified with the antineutrino iV£ = (1, 1, 0). It carries lepton number 
equal to L = — 1. It is possible to introduce a new scalar field 5 4 which may couple to N1u c L 
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and d^d^, thus violating the baryon number. It is remarkable that the choice for the lepton 
number of iV£ leads to no new gauge boson which violates B — L. These considerations 
do not apply to supersymmetric models though (for a review see ]5q| ). Indeed, for every 
fermionic degree of freedom there exist a superpartner (squark or slepton) which does have 
the same quantum number. Furthermore, in the Minimal Supersymmetric Standard Model 
(MSSM) one has to introduce two Higgs doublets Ht = (1,2,-1/2) and H 2 = (1,2,1/2) 
and the corresponding fermionic superpartners, the so-called higgsinos Hi 2- Finally, every 
gauge boson has its own superparner, the gaugino. In this large zoo of new particles, one 
can easily find couplings that violate B and B — L. For instance, the higgsino Hi may couple 

10 — ft 

to the quark doublet Ql and to the scalars S and S of the Table 2. The pair H[Q' L has 
baryon number B = —1/3 and B — L = —1/3 and both quantum numbers are not conserved. 
Nevertheless, in the supersymmetric models which are phenomenologically acceptable, even 
without considering the presence of superheavy particles, it is necessary to suppress some 
supersymmetric couplings which would lead at the weak scale to a proton decay at a rate 
which is too fast for being in agreement with the tight experimental constraints . One 
commonly accepted solution is to introduce a discrete symmetry Z%, called -R-parity, under 
which all the fields of the SM are even and all the superpartners are odd. The scalar 
component of any chiral supermultiplet has the following i?-parity number 

R=(-lf( B - L \ (101) 

while the corresponding fermion has the same number multiplied by —1. If we impose that 
/2-parity is exact, then it is easy to check that, besides suppressing the fast proton decay 
at the weak scale, one avoids the presence B and L violating couplings of heavy fields with 
the light fermionic fields of the MSSM. Indeed, all the heavy bosons of Tables 1 and 2 
have i?-parity R = 1, while - for instance - the pair H\Q^ L has parity R = — 1. Similar 
considerations apply to other fermionic pairs. 

We conclude that in the GUTs, both supersymmetric and non-supersymmetric, no B — L 
asymmetry may be generated through the out-of-equilibrium decay of gauge boson fields. 
We will mention in the following - though - that the generation of such an asymmetry is 
possible in the framework of (supersymmetric) SO(W) via the out-of-equilibrium decay of 
the right-handed (s)neutrino, i.e. via the decay of a superheavy fermion (scalar). 

After having learned that GUTs are the perfect arena for baryon number violating 
interactions, we will illustrate now some features of the out-of-equilibrium decay scenario 
within some specific GUTs, like SU(5) and SO(10). 



4.3.1 The case of SU(5) 

The gauge group SU(5) is the smallest group containing the SM gauge group and as such 
it represents the most appealing candidate to build up a Grand Unified Theory. The non- 
supersymmetric version of SU (5) is - however - already ruled out by its prediction of 
the proton lifetime r v ~ 10 30 years, which is in disagreement with the experimental lower 
bound t p > 10 32 years @. Recent precise measurements of coupling constants at LEP 
suggest that the supersymmetric extension of SU(5) gives a consistent picture of coupling 
unification |^, 39, and is a viable possibility. 
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The fermionic content of SU(5) is the same as the one in the SM. Therefore, as we 
explained in Section 4.3, it is not possible to create any asymmetry in B — L. Fermions are 
assigned to the reducible representation 5/ © 10/ as 

5 f = {d c L ,£ L } (102) 

and 

10/ = {Q L ,ul,e c L } . (103) 

There are 24 gauge bosons which belong to the adjoint representation 24y and may couple 
to the fermions through the couplings 

^=2A V [(5/)t (5/) + (10/)t 10/] . (104) 

Among the 24 gauge bosons there are the bosons XY = V 1 = (3, 2, —5/6) (and their CP- 
conjugate) which may decay violating the baryon number: XY — > QL, QQ, where Q and L 
denotes an arbitrary quark and lepton, respectively. They have electric charges Qx = — 1/3 
and Qy = —4/3. The mass and the couplings of these bosons are determined by the gauge 
coupling unification 

Mxy — 5 x 10 1 GeV, acuT — 1/45, non — supersymmetric SU(5), 

M X y ^ 10 16 GeV, ogut^I/24, supersymmetric SU(5). (105) 

While in the gauge sector the structure is uniquely determined by the gauge group, in the 
Higgs sector the results depend upon the choice of the representation. The Higgs fields 
which couple to the fermions may be in the representation 5// or in the representations 
10//, 15//, 45// and 50//. If we consider the minimal choice 5//, we obtain 

ku (10/) T (10/) 5// + h D (5/) T (10/) 5//, (106) 

where hu,D are matrices in the flavor space. The representation 5// contains the Higgs 
doublet of the SM, (1,2,1/2) and the triplet S 1 = (3,1,-1/3) which is B-violating. Un- 
fortunately, this minimal choice of the Higgs sector does not suffice to explain the baryon 
number of the Universe. The CP violation is due to the complex phases which cannot be 
reabsorbed by field redifinition (they are physical) in the Yukawa sector. At the tree-level 
these phases do not give any contribution to the baryon asymmetry and at the one-loop 
level the asymmetry is proportional to 

ImTr (h\jhuh\jh D ) =0, (107) 

where the trace is over generation indices. This is because the Higgs on the external and 
internal legs of the one-loop interference diagrams is the same. A net baryon number only 
appears at three-loop, resulting in a baryon asymmetry ~ 10~ 16 which is far too small to 
explain the observed one. The same problem in present in the supersymmetric version of 
SU (5) where one has to introduce two Higgs supeerfields 5// and 5jj . 

The problem of too tiny CP violation in SU(5) may be solved by complicating further 
the Higgs sector. One may introduce an extra scalar 5' H with the same quantum numbers 
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of 5h, but with a different mass and/or lifetime [101 1 . In that case one- loop diagrams with 
exchange of 5' H instead of 5# can give rise to a net baryon number proportional 



ImTr (h'^huh'^hD 



(108) 



where h' v D are the couplings of 5'^ to QL and QQ, respectively. A second alternative is 
to introduce a different second Higgs representation. For example, adding a Higgs in the 
45 representation of SU(5) an adequate baryon asymmetry may be producedb for a wide 
range of the parameters |oT 



4.3.2 The case of 50(10) 

In the GUT based on 50(10) the spontaneous breaking down to the SM gauge group is 



generally obtained through different steps (for a general review, see [117])- The main two 
channels are 



5O(10) 
5O(10) 



SU(3)c®SU(2) L ®U(l] 



y ■> 



Mqut n Mr „ M c ^ M B _ L 
— ► Cx224 — > 1*214 —> t*2113 ~ ► 

M ^ UT G 2 24 -° G2213 ^ Gam L SU(3) C ® SU(2) L ® U(l) Y , (109) 



where 



G 22 4 = SU(2) L ® SU(2) R ®SU{A), 

G214 = SU{2) L ® U{l) hR ® 5C/(4), 

G2113 = SU(2) L ®U(l) hR ®U(l) B - L ®SU(3) c , 

G2213 = SU(2) L SU(2) R ® tf(l) fl -L ® 5Z7(3)c, 



(110) 



where the four intermediate scales have not to be necessarily different from each other. We 
notice that, 

- if we are interested in the generation of an asymmetry in B — L, the relevant scale is 
the scale at which the abelian group U(1)b-l breaks down, i.e. Mb-l, arid not the Grand 
Unification scale Mqut; 

- it is usually not possible to generate any baryon asymmetry at the scale Mqut- 
Indeed, the fermionic content of 5O(10) is the one of the SM plus a right-handed neutrino 
iV£ = (1,1,0). All the fermions belonging to the same generation are contained in the 
spinorial representation 16 1. Differently from what happens for the case of SU(5), now all 
the fermions posses the corresponding antifermions and it is possible to define a conjugation 
operator of the charge C starting from the operators of 5O(10), in such a way that, if 
5O(10) is not broken, then C is conserved |9(|. In the simplest mechanism for the breaking 
of 5O(10), the one into G224, a crucial role is played by the 54#. In such a case there is 
a symmetry SU(2)l 



SU(2)r [|104| , |93| , |94fl with equal coupling constants gi and gn and 
consequently, C is still a symmetry of the theory. It is possible to see that, with this choice 
of the Higgs representation, C is not broken until U(1)b-l is broken, i.e. at the scale Mb-l 
0. At this scale, the right-handed neutrino acquires a Majorana mass M/v = 0(M b ~l) and 

3 In fact, if one uses the 210.tr to break 50(10), one can maintain the gauge part of the left-right symmetry 
but not C. Th eimplications for baryogenesis are discussed in p3j|. 
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its out-of-equilibrium decays may generate a nonvanishing B — L asymmetry . We will 
return to this point later. With a more complicated choice of the Higgs representation it is 
possible to break C at the scale Mr where SU (2)# is broken and in such a case baryogenesis 
may take place at that scale. 

5 The out-of-equilibrium decay scenario and the thermal his- 
tory of the Universe 

The out-of-equilibrium scenario that we have depicted in the previous section is operative 
only if a nonequilibrium number density of X heavy bosons was present in the early Universe. 
Usually massive particles are in equilibrium at at high temperatures, T ^> Mx and their 
number density exceeds the equilibrium one when T becomes of the same order of the mass 
Mx- We have seen that, if the decay rate is small enough around T ~ Mx, see Eq. (J58l), 
then departure from equilibrium is attained and the subsequent decays of X and X particles 
may produce the observed baryon number asymmetry. The basic assumption - however - 
of this picture is that the superheavy bosons were as abundant as photons at very high 
temperatures T > Mx- 

If the X particles are gauge or Higgs bosons of Grand Unification, the situation is 
somewhat more complicated because they might have never been in thermal equilibrium 
at the very early stages of the evolution of the Universe. Even if the temperature of the 
primeval plasma was higher than the Grand Unified scale Mgut ~ 10 16 GeV, the rate of 
production of superheavy particles would be smaller than the expansion rate of the Universe 
and the number density of superheavy bosons could always be smaller than the equilibrium 
one. Secondly, the temperature of the Universe might be always smaller than Mgut and 
correspondingly the thermally produced X bosons might be never as abundant as photons, 
making their role in baryogensis negligible. All these considerations depend crucially upon 
the thermal history of the Universe and deserve a closer look. 



5.1 Inflation and reheating: the old days 

The flatness and the horizon problems of the standard big bang cosmology are elegantly 
solved if during the evolution of the early Universe the energy density happened to be 
dominated by some form of vacuum energy and comoving scales grow quasi-exponentially 



[54]. An inflationary stage is also required to dilute any undesirable topological defects left 
as remnants after some phase transition taking place at early epochs. 

The vacuum energy driving inflation is generally assumed to be associated to the poten- 
tial V(4>) of some scalar field (ft, the inflaton, which is initially displaced from the minimum 
of its potential. As a by-product, quantum fluctuations of the inflaton field may be the 
seeds for the generation of structure and the fluctuations observed in the cosmic microwave 
background radiation, 5T/T ~ 10 -5 M, M M. 



Inflation ended when the potential energy associated with the inflaton field became 
smaller than the kinetic energy of the field. By that time, any pre-inflation entropy in the 
Universe had been inflated away, and the energy of the universe was entirely in the form of 
coherent oscillations of the inflaton condensate around the minimum of its potential. The 
Universe may be said to be frozen after the end of inflation. We know that somehow the 
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low-entropy cold Universe dominated by the energy of coherent motion of the eft field must 
be transformed into a high-entropy hot Universe dominated by radiation. The process by 
which the energy of the inflaton field is transferred from the inflaton field to radiation has 
been dubbed reheating. In the old theory of reheating p5| , [p, the simplest way to envision 
this process is if the comoving energy density in the zero mode of the inflaton decays into 
normal particles, which then scatter and thermalize to form a thermal background. It is 
usually assumed that the decay width of this process is the same as the decay width of a 
free inflaton field. 

Of particular interest is a quantity known as the reheat temperature, denoted as Trh- 
The reheat temperature is calculated by assuming an instantaneous conversion of the energy 
density in the inflaton field into radiation when the decay width of the inflaton energy T^, 
is equal to H , the expansion rate of the universe. 

The reheat temperature is calculated quite easily. After inflation the inflaton field ex- 
ecutes coherent oscillations about the minimum of the potential. Averaged over several 
oscillations, the coherent oscillation energy density redshifts as matter: p^ oc a -3 , where 
a is the Robertson-Walker scale factor. If we denote as pi and aj the total inflaton en- 
ergy density and the scale factor at the initiation of coherent oscillations, then the Hubble 
expansion rate as a function of a is 

Equating H{a) and leads to an expression for aj/a. Now if we assume that all available 
coherent energy density is instantaneously converted into radiation at this value of ai/a, 
we can find the reheat temperature by setting the coherent energy density, p^ = pj(ai/a) 3 , 
equal to the radiation energy density, pr = (tt 2 / '30) g^T^ H , where <?* is the effective number 
of relativistic degrees of freedom at temperature Trh- The result is 



r - = (s^:) 1/4 ^= ' 2 (™) l/4 ^- 

In the simplest chaotic inflation model, the inflaton potential is given by 



:n2) 



1 



V(cP) = ^M|</> 2 , (113) 



2 

with 

~ 10 13 GeV (114) 

in order to reproduce the observed temperature anisotropics in the microwave background 
[M. Writing = a^M^, one finds 

T^~10 15 ^GeV. (115) 

5.2 GUT baryogenesis and the old theory of reheating: a Herculean task 

There are very good reasons to suspect that GUT baryogenesis is not in a good shape in 
the old theory of reheating. 
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5.2.1 Kinematical suppression of superheavy particles 

The density and temperature fluctuations observed in the present universe, 5T/T ~ 10 -5 , 
require the inflaton potential to be extremely flat - that is < 1. This means that the 
couplings of the inflaton field to the other degrees of freedom cannot be too large, since 
large couplings would induce large loop corrections to the inflaton potential, spoiling its 
flatness. As a result, Trh is expected to be much smaller than 10 14 GeV by several orders 
of magnitude. As we have seen, the unification scale is generally assumed to be around 10 16 
GeV, and £>-violating gauge bosons should have masses comparable to this scale. Baryon- 
number violating Higgs bosons may have a mass one or two orders of magnitude less. For 
example, in SU (5) the B violating Higgs bosons in the five-dimensional representation that 
may have a mass as small as 10 14 GeV. In fact, these Higgs bosons are more likely than 
gauge bosons to produce a baryon asymmetry since it is easier to arrange the requisite CP 
violation in the Higgs decay. But even the light U-violating Higgs bosons are expected 
to have masses larger than the inflaton mass, and it would be kinematically impossible to 
create them directly in <f> decay, (j) — ► XX. This is because one expects 

M «M X . (116) 



5.2.2 Thermal production of heavy particles 

One might think that the X bosons could be created by thermal scattering during the stage 
of thermalization of the decay products of the inflaton field. Indeed, the reheat temperature 
is best regarded as the temperature below which the Universe becomes radiation dominated. 
In this regard it has a limited meaning. For instance, it should not be interpretated as the 
maximum temperature obtained by the universe during reheating. The maximum temper- 
ature is, in fact, much larger than Trh- One implication of this is that it is incorrect, 
to assume that the maximum abundance of a massive particle species X produced after 
inflation is suppressed by a factor of exp(— Mx/Trh) lH] and therefore it is incorrect to 
conclude that GUT baryogenesis is imcompatible with models of inflation where the re- 
heating temperature is much smaller than the GUT scale and, in general, than the mass of 
the X particles, Trh <C Mx- Particles of mass much greater than the eventual reheating 
temperature Trh may be created by the thermalized decay products of the inflaton. In- 
deed, a stable particle species X of mass Mx would be produced in the reheating process in 
sufficient abundance that its contribution to closure density today would be approximately 
M x (a\v\) (g* / 10 3 ) (Mx /W 4 Trh) , where is the number of effective degrees of freedom of 
the radiation energy density and (cr|u|) is the thermal average of the X annihilation cross 
section times the M0ller flux factor. Thus, particles of mass as large as 10 4 times the re- 
heating temperature may be produced in interesting abundance [26|. The number density 
nx of particles X after freeze out and reheating may be easily computed [26| 

n 7 \gj \Mx) \M X ) K ' 
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and is not exponentially suppressed. It is easy to check that for such small values of Tru, 



the ratio (117) is always much larger than the equilibrium value 




= («3)j(2fwJ (U8> 

This result is crucial for the out-of-equilibrium decay scenarios of baryogenesis. For instance, 
as we shall see, in theories where B — L is a spontaneously broken local symmetry, as 
suggested by 50(10) unification, the cosmological baryon asymmetry can be generated 
by the out-of-equilibrium decay of the lightest heavy Major ana right-handed neutrino iVf , 
whose typical mass is about 10 10 GeV [fj^| . For reheat temperatures of the order of 10 9 
GeV, the number density of the right-handed neutrino is about 3 x 10~ 2 n 7 and one can 
estimate the final bayon number to be of the order of B ~ (n^/n^)(e/ 'g*) — 10~ 4 e, where 
e is the coefficient containing one-loop suppression factor and CP violating phases. The 
observed value of the baryon asymmetry, B ~ 10 -10 , is then obtained without any fine 
tuning of parameters. 

Exercise 2 

Compute the maximum temperature during the process of reheating. Hint: Consider 
the early-time solution for radiation (i.e. when H 3> r<j, and before a significant fraction of 
the comoving coherent energy density is converted to radation). 

5.2.3 The gravitino problem 

There is one more problem associated with GUT baryogenesis in the old theory of reheat- 
ing, namely the problem of relic gravitinos [41]. If one has to invoke supersymmetry to 
preserve the flatness of the inflaton potential, it is mandatory to consider the cosmologi- 
cal implications of the gravitino - a spin-(3/2) particle which appears in the extension of 
global supersymmetry to local supersymmetry - or supergravity [pC| J. The gravitino is the 
fermionic superpartner of the graviton and has interaction strength with the observable 
sector - that is the SM particles and their superpartners - inversely proportional to the 
Planck mass. One usually associates the scale of supersymmetry breaking with the elec- 
troweak scale in order to handle the hierarchy problem [103] and the mass of the gravitino 



is of order of the weak scale, 7713/2 = 0(1) TeV. The decay rate of the gravitino is given by 

r^t^^&tt- (119> 

The slow decay rate of the gravitinos is the essential source of the cosmological problems 
because the decay products of the gravitino will destroy the 4 He and D nuclei by photodis- 
sociation, and thus successful nucleosynthesis predictions. The most stringent bound comes 
from the resulting overproduction of D + 3 He, which would require that the gravitino abun- 
dance is smaller than ~ 10~ 10 relative to the entropy density at the time of reheating after 
inflation p0[ 

— & (10~ 10 - HT 11 ). (120) 
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The Boltzmann equation governing the number density of gravitinos n 3 / 2 during the ther- 
malization stage after inflation is 



dn 3/2 
dt 



+ 3-ffn 3/2 ~ (Etot^njjght, (121) 



where S tot oc 1/Mp is the total cross section determining the rate of production of gravitinos 
and riiight ~ T 3 represents the number density of light particles in the thermal bath. The 



number density of gravitinos at thermalization is readily obtained solving Eq. ( 121 ) and 
reads 

^Il^ l0 -2l3H. (122) 

s M P y ' 



Comparing Eqs. ( |120| ) and (122), one may obtain an upper bound on the reheating tem- 
perature after inflation 

Trh < (10 10 - 10 11 ) GeV. (123) 

Therefore, if Trh ~ -Mgut, gravitinos would be abundant during nucleosynthesis and 
destroy the good agreement of the theory with observations. However, if the initial state 
after inflation was free from gravitinos, the reheating temperature seems to be too low to 
create superheavy X bosons that eventually decay and produce the baryon asymmetry - 
even taking into account the previous considerations about the fact that the maximum 
temperature during reheating is not Trh [[^ |26| . 

5.3 Inflation and reheating: the new wisdom 

The outlook for GUT baryogenesis has brightened recently with the realization that reheat- 



ing may differ significantly from the simple picture described above |68|, S3, |34|, |6^, p6|. In 
the first stage of reheating, called preheating J6S|], nonlinear quantum effects may lead to 
an extremely effective dissipational dynamics and explosive particle production even when 
single particle decay is kinematically forbidden. Particles can be produced in the regime 
of a broad parametric resonance, and it is possible that a significant fraction of the energy 
stored in the form of coherent inflaton oscillations at the end of inflation is released after 
only a dozen or so oscillation periods of the inflaton. What is most relevant for these lec- 



tures is that preheating may play an extremely important role for baryogenesis [74, H, 75] 



and, in particular, for GUT generation of the baryon asymmetry. Indeed, it was shown in 



[74, 75 1 that the baryon asymmetry can be produced efficiently just after the preheating era, 
thus solving many of the problems that GUT baryogenesis had to face in the old picture of 
reheating. 

The presence of a preheating stage at the beginning of the reheating process is based 
on the fact that, for some parameter ranges, there is a new decay channel that is non- 
perturbative: due to the coherent oscillations of the inflaton field stimulated emissions 



of bosonic particles into energy bands with large occupancy numbers are induced [68]. 
The modes in these bands can be understood as Bose condensates, and they behave like 
classical waves. The back-reaction of these modes on the homogeneous inflaton field and 
the rescattering among themselves produce a state that is far from thermal equilibrium and 
may induce very interesting phenomena, such as non-thermal phase transitions |]69| , 121 , |112| 



with production of a stochastic background of gravitational waves [66] and of heavy particles 
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in a state far from equilibrium, which may constitute today the dark matter in our Universe 



The idea of preheating is relatively simple, the oscillations of the inflaton field induce 
mixing of positive and negative frequencies in the quantum state of the field it couples to 
because of the time- dependent mass of the quantum field. Let us focus - for sake of simplicity 
- to the case of chaotic inflation, with a massive inflaton cj> with quadratic potential V(4>) = 
^M 2 (f) 2 , Mff, ~ 10 13 GeV, and coupled to a massless scalar field x y i a the quartic coupling 

The evolution equation for the Fourier modes of the \ field with momentum k is 

X k + to 2 k X k = 0, (124) 

with 

X k = a 3 / 2 (t) Xk , 

col = fc 2 /a 2 W + <?V(i). (125) 
This Klein-Gordon equation may be cast in the form of a Mathieu equation 

X'l + [A{k) - 2q cos 2z]X k = 0, (126) 

where z = M^t and 

2 



" = ^ (127 > 

where $ is the amplitude and is the frequency of inflaton oscillations, <p(t) = $(t) sm(M^,t). 
Notice that, at least initially when $ = cMp < Mp 

gl^L „ g 2 c 2fM| „ g 2 c 2 x 1Q 12 ^ 1 (128) 

M 4M| M? v ' 

and the resonance is broad. For certain values of the parameters (^4, q) there are exact 
solutions X k and the corresponding number density n k that grow exponentially with time 
because they belong to an instability band of the Mathieu equation (for a recent compre- 
hensive review on preheating after chaotic inflation [f7C| l and references therein) 

X k oc e^ kM ^ ^n k (x e 2WsM ^, (129) 

where the parameter fi k depends upon the instability band and, in the broad resonance 
case, q 3> 1, it is ~ 0.2. 

These instabilities can be interpreted as coherent "particle" production with large oc- 
cupancy numbers. One way of understanding this phenomenon is to consider the energy of 
these modes as that of a harmonic oscillator, E k = \X k \ 2 /2 + uj^\X k \ 2 /2 = 0J k {n k + 1/2). 
The occupancy number of level k can grow exponentially fast, n k ~ exp(2/ifcM^t) 1, 
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and these modes soon behave like classical waves. The parameter q during preheating de- 
termines the strength of the resonance. It is possible that the model parameters are such 
that parametric resonance does not occur, and then the usual perturbative approach would 
follow, with decay rate T^. In fact, as the Universe expands, the growth of the scale factor 
and the decrease of the amplitude of inflaton oscillations shifts the values of (A, q) along the 
stability /instability chart of the Mathieu equation, going from broad resonance, for g > 1, 
to narrow resonance, q <C 1, and finally to the perturbative decay of the inflaton. 

It is important to notice that, after the short period of preheating, the Universe is 
likely to enter a long period of matter domination where the biggest contribution to the 
energy density of the Universe is provided by the residual small amplitude oscillations of 
the classical inflaton field and /or by the inflaton quanta produced during the back-reaction 
processes. This period will end when the age of the Universe becomes of the order of 
the perturbative lifetime of the inflaton field, t ~ 177 . At this point, the Universe will 
be reheated up to a temperature Trh given in ( |112|) obtained applying the old theory of 
reheating described in the previous section. 



5.4 GUT baryogenesis and preheating 

A crucial observation for baryogenesis is that even particles with mass larger than that of 
the inflaton may be produced during preheating. To see how this might work, let us assume 
that the interaction term between the superheavy bosons and the inflaton field is of the 
type g 2 cf) 2 \X\ 2 . During preheating, quantum fluctuations of the X field with momentum k 
approximately obey the Mathieu equation where now 

k 2 + M 2 

A(k) = k -±^ + 2q. (130) 

Particle production occurs above the line A = 2q. The width of the instability strip scales 
as q^l 2 for large q, independent of the X mass. The condition for broad resonance p8| , fF3| 

A-2q<q l l 2 (131) 

becomes 

k 2 + Ml $ . > 

— o < 9 , 132) 

which yields for the typical energy of X bosons produced in preheating 

E 2 x = k 2 + M] c <g^M 4> , (133) 

By the time the resonance develops to the full strength, <I> 2 ~ 10 -5 Aip. The resulting 
estimate for the typical energy of particles at the end of the broad resonance regime for 
M ~ 10 13 GeV is 

E x ~ 10-y /2 y / AfyMp ~ <? 1/2 10 15 GeV. (134) 

Supermassive X bosons can be produced by the broad parametric resonance for Ex > Mx, 
which leads to the estimate that X production will be possible if Mx < g 1//2 10 15 GeV. 
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For g 2 ~ 1 one would have copious production of X particles (in this regime the problem 
is non-linear from the beginning and therefore g 2 = 1 has to be understood as a rough 
estimate of the limiting case) as heavy as 10 15 GeV, i.e., 100 times greater than the inflaton 
mass. The only problem here is that for large coupling g, radiative corrections to the effective 
potential of the inflaton field may modify its shape at cj> ~ Mp\. However, this problem 
does not appear if the flatness of the inflaton potential is protected by supersymmetry. 

This is a significant departure from the old constraints of reheating. Production of X 
bosons in the old reheating picture was kinematically forbidden if < Mj, while in the 
new scenario it is possible because of coherent effects. It is also important to note that the 
particles are produced out-of-equilibrium, thus satisfying one of the basic requirements to 
produce the baryon asymmetry [ 114 ]. 

Scattering of X fluctuations off the zero mode of the inflaton field limits the maxi- 
mum magnitude of X fluctuations to be (X 2 ) max m M^/g 2 [p5| . For example, (X 2 ) max ~ 
10~ 10 M| in the case Mx = 10 Ma,. This restricts the corresponding number density of 
created X-particles. 

A potentially important dynamical effect is that the parametric resonance is efficient only 
if the self-interaction couplings of the superheavy particles are not too large. Indeed, a self- 
interaction term of the type A] A] 4 provides a non-thermal mass to the X boson of the order 
of (X(X 2 )) 1 ^ 2 , but this contribution is smaller than the bare mass Mx, if A < g 2 M x /M 2 . 
Self-interactions may also terminate the resonance effect because scattering induced by the 
coupling A may remove particles from the resonance shells and redistribute their momenta 
p8| . But this only happens if, again, A >j 2 [64]. 

The parametric resonance is also rendered less efficient when the X particles have a 
(large) decay width Tx, which is essential for the out-of-equilibrium decay to take place. 
Roughly speaking, one expects that the explosive production of particles takes place only 
if the typical time, r e , during which the number of X bosons grows by a factor of e, is 
smaller than the decay lifetime tx = . During the broad resonance regime, typically 
T e < 10 M^ 1 . If we write the decay width by Tx = olx Mx, this requires ax < O.IM^/Mx- 
Notice that smaller values of Tx are favored not only because particle production is made 
easier, but also because the superheavy particles may remain out-of-equilibrium for longer 
times, thus enhancing the final baryon asymmetry. 

Using the methods developed in Refs. [63, 64, 65], one can study numerically the pro- 
duction of massive, unstable X particles in the process of the inflation decay [fTq] ■ Let us 
consider a model in which the oscillating inflaton field cj) interacts with a scalar field X 
whose decays violate baryon number B. As we have learned, the simplest possibility for the 
X-particle is the Higgs field in the five-dimensional representation of SU(5). We assume 
standard kinetic terms, minimal coupling with gravity, and a very simple potential for the 
fields of the form 

V(<f>, X) = hl 2 <p 2 + \m\X 2 + \g 2 ^X 2 . (135) 

A fundamental parameter in GUT baryogenesis is nx, the number density of the super- 
massive leptoquarks whose decays produce the baryon asymmetry. It will depend upon the 
value of r and q. 

Since the supermassive bosons are more massive than the inflaton, one expects small 
kinetic energy in the excitations of the X field. From the potential of Eq. (|135|), the square 
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Figure 4: The variance of X with model parameters q = 10 6 , m x = 2, and T = Tx/M^ = 
6 x 10~ 2 is shown by the lower solid curve as a function of time. The upper solid curve 
corresponds to the inflaton zero mode. The dotted curves represent the same quantities for 
r = 0. 



of the effective mass of the X field is 

rEFF\2 



(M 



X 



M 2 x +g 2 {<f> 2 ) 



and the energy density in the X field will be 

Px^(M 2 x +g 2 {<t ) 2 )){X 2 ). 



(136) 
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Writing 



as 



+ 



), one can define an analog of the X-particle number density as 



n x = Px/M x 



EFF 



Aq(4>l + <^ 2 »/0 2 (O) + m 2 ] V2 M^X 2 ), 



(138) 



where m x = Mx/M^. 



Eq. (|13q ) enables one to calculate the number density of the created X-particles if the 
variances of the fields, (X 2 ), (5(f) 2 ), and the inflaton zero mode ^o( r ) (here t and r are 
related by M^dt = a{T)dr) are known. 

The time evolution of the variance, (X 2 ), and of the inflaton zero mode, (cp), is shown 
in Fig. 4, by the solid curves for the case q = 10 6 , m x = 2, and T = 6 x 10~ 2 . We see that 
the particle creation reaches a maximum at r « 10.8 when (X 2 ) ps 1CT 9 in the "valleys" 
between the peaks. At later times, r > 10.8, particle creation by the oscillating inflaton 
field can no longer compete with X-decays due to the non-zero value of T. For comparison, 



we show in the same figure the case T = represented by the dotted curves [35|. In the 
r = case, particle creation is able to compete with the expansion of the universe so that 
(X 2 ) remains roughly constant. 
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Figure 5: The maximum value of the variance of the X-field, (X^), is shown as function 
of r. Stars mark (X^) obtained in the full non-linear problem. {X 2 nax ) in the Hartree 
approximation is shown by the dotted curve for q = 10 4 , m x = 0.1, and by the solid curve 
for q = 10 6 , m x = 2. 



Using Eq. ( |138[) , one finds for the maximum number density of created X-particles 



n x 



4#§(10.8)/^j (0) + ml] 1/2 M^X 



10 3 + m 2 x 



1/2 



M^X 2 ) «30M^(X 2 ). (139) 



It is easy to understand that if we increase the value of T, the parametric resonance will 
not be able to compete with the decay of X at earlier times. Moreover, for sufficiently large 
values of T, the resonance will be shut off in the linear regime. 

In exploration of parameter space it turns out more convenient to go to the Hartree 
approximation which requires much less computing resources. The maximum value of the 
variance of X reached during the time evolution of the fields in the Hartree approximation 
is shown in Fig. || as a function of the parameters of the model. Here the stars also show 
the maximum of (X 2 (t)) in the full non-linear problem for a few values of T. At small T 
the Hartree approximation overestimates (X 2 ) significantly |63, |65]] . Nonetheless, at large 
values of T it is a quite reliable approach. One may see that (X 2 ) drops sharply when 
T > 0.2, and this critical value of T does not depend significantly upon mx or q Iffq] , 

The most relevant case with q = 10 s , where X-bosons as massive as ten times the 
inflaton mass can be created, is shown in Fig. 6 in the Hartree approximation. Note, that 
two lower curves which correspond to T equal to 0.08 and 0.12 never reach the limiting value 
(X 2 ) max ~ 10~ 10 Mp, which is imposed by rescattering [65], and the Hartree approximation 
ought to be reliable in this cases. 

As outlined above, one may consider a three part reheating process, with initial con- 
ditions corresponding to the frozen universe at the end of inflation. The first stage is 
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Figure 6: The time dependence of the variance of X in the Hartree approximation with 



model parameters q 
0.08, 0.12. 



10 s , 



10 and for three values of T, from top to bottom: 0.04, 



explosive particle production, where a fraction 5 of the energy density at the end of pre- 
heating is transferred to X bosons, with (1 — 5) of the initial energy remaining in cj> coherent 
oscillation energy. We assume that this stage occurs within a few Hubble times of the end 
of inflation. The second stage is the X decay and subsequent thermalization of the decay 
products. We assume that decay of an X—X pair produces a net baryon number e, as well 
as entropy. Reheating is brought to a close in the third phase when the remaining energy 
density in <j> oscillations is transferred to radiation. 

The final baryon asymmetry depends linearly upon the ratio 5 between the energy stored 
in the X particles at the end of the preheating stage and the energy stored in the inflaton 
field at the beginning of the preheating era [74]. The description simplifies if we assume 
zero initial kinetic energy of the Xs. One may also a assume that there are fast interactions 
that thermalize the massless decay products of the X. Then in a co- moving volume a 3 , the 



total number of X bosons, Nx = nxa , the total baryon number, Nb 
dimensionless radiation energy, R = prcl 4 , evolve according to 



ubcl , and the 



N x 
N B 



-Y x (n x - N* Q 



R = -aM x N x ; 



-eNx-TxNslN^/No 



(140) 



N^Q is the total number of Xs in thermal equilibrium at temperature T oc i? 1//4 , and iVo is 
the equilibrium number of a massless degree of freedom in a comoving volume. 



37 



-I — I I 1 1 1 1 1 1 1 — I I 1 1 1 1 1 1 1 — I I 1 1 1 1 1 1 1 — I I I 1 1 1 1 1 1 — I I 1 1 1 1 1 1 1 — I I 1 1 1 1 1 1 1 — I I I 1 1 1 1 1 1 — I I 1 1 

•-, (P.a- 3 )/^' ar») 1 




B/el 



(Pxa-VWaf 3 ) 



^3/2 ~i 



J ■ c I 



10 1 10 2 10 3 10 4 10 5 10 8 10 7 10 8 

a/a, 



Figure 7: The evolution of the baryon number, the X number density, the energy density 
in (j> oscillations, and the gravitino-to-entropy ratio as a function of the scale factor a. 



Fig. 7 shows the results of an integration of Eqs. ( |140| ) in a toy model with M, 



10 ld GeV, M x =_ 10 i4 GeV, T x = 5 x KT b M x , ^ = 5 x 10 -iU M^ , and two degrees of 
freedom (b and b). Initial conditions were chosen at a = ai to be px = Pcj> ~ 10~ 4 M|Mp, 
and R = Nb = 0. The px = P<j> assumption corresponds to 5 = 1/2. The baryon number 
B = ub/s rapidly rises. However B decreases as entropy is created and X inverse reactions 
damp the baryon asymmetry. After most of the energy is extracted from the initial X 
background, the baryon number is further damped as entropy is created during the decay 
of energy in the cj> background. One can also numerically integrate the equation governing 
the number density of gravitinos ra 3/ / 2 - The result for G 3 / 2 = n^^/s is shown in Fig. 7. 
Notice that, even though gravitinos are copiously produced at early stages by scatterings 
of the decay products of the X, G 3 / 2 decreases as entropy is created during the subsequent 
decay of energy in the cfi background. 

Since the number of X bosons produced is proportional to 5, the final asymmetry is 
proportional to 5 and B/e ~ 10 -9 can be obtained for 5 as small as 10 -6 . One can estimate 
this ratio as 

6 fg-_ (X^ 



3 x 10 



106 ™ x Ml 



(141) 



Therefore, for q = 10 8 and m x = 10, 5 is of the order of 3 x 10 8 (X 2 )/Mp. Since the final 
baryon asymmetry scales approximately as T _1 and is given by 5 ~ SxlO^etr/SxlO" 5 )- 1 
|f74}| , where e is an overall parameter accounting for CP violation, one can see that the 
observed baryon asymmetry B~4x 10~ n may be explained by the phenomenon of GUT 
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baryogenesis after preheating if 



M| \ e j V5 x IO- 5 

From Fig. 6 we can read that this only may happen if 

T x < 1(T 3 M x . (143) 



This result may be considered very comfortable since we can conclude that whenever the res- 
onance develops, i.e., when T x < lO^M^ = KT 2 M X , GUT baryo genesis after preheating 
is so efficient that the right amount of baryon asymmetry is produced for almost the entire 
range of values of the decay rate Tx- In other words, provided that superheavy X-bosons 
are produced during the preheating stage, they will be ineffective in producing the baryon 
asymmetry only if their decay rate falls in the range 10~ 3 Mx < Fx < 10~ 2 Mx- GUT 
baryogenesis after preheating solves many of the serious drawbacks of GUT baryogenesis in 
the old theory of reheating where the production of superheavy states after inflation was 
kinematically impossible. Moreover, the out-of-equilibrium condition is naturally attained 
in our scenario since the distribution function of the X-quanta generated at the resonance 
is far from a thermal distribution. This situation is considerably different from the one 
present in the GUT thermal scenario where superheavy particles usually decouple from the 
thermal bath when still relativistic and then decay producing the baryon asymmetry. It is 
quite intriguing that out of all possible ways the parametric resonance may develop, Nature 
might have chosen only those ways without instantaneous thermalization and also with a 
successful baryogenesis scenario. 



6 The baryon number violation in the Standard Model 

In this section we will be concerned with the violation of the baryon and lepton number in 
the SM. It is well-known that by considering the most general Lagrangian invariant under 
the SM gauge group SU(3)c <8> SU(2)l (g> U(l)y and assuming that the Higgs fields are 
color singlets, the Lagrangian is automatically invariant under global abelian symmetries 
which may be identified with the baryonic and leptonic symmetries. They are - therefore - 
accidental symmetries. As a result, it is not possible to violate B and L at the tree- level and 
at any order of perturbation theory: the proton is stable in the SM and any perturbative 
process which violates B and/or L in Grand Unified Theories is necessarily suppressed by 
powers of Mqut/Mw- Nevertheless, in many cases the perturbative expansion does not 



describe all the dynamics of the theory and - indeed - in 1976 't Hooft [120] realized that 
nonperturbative effects (instantons) may give rise to processes which violate the combination 
B + L, but not the orthogonal combination B — L. The probability of these processes to 
occur is exponentially suppressed, ~ exp(— Air /aw) ~ 10~ 150 where aw = ffl/^TT is the 
weak gauge coupling, and probably irrelevant today. In more extreme situations - like the 



primordial Universe at very high temperatures [|3J, 67, [76J - baryon and lepton number 
violation processes may be fast enough to play a significant role in baryogenesis. This will 
be the subject of the present section. 
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6.1 The B + L anomaly 



The violation of the baryonic number within the SM is due to the fact that the current 
corresponding to the global abelian group Ub+l~ even though it is conserved at the classical 
level - is not conserved at the quantum level, that is the Ub+ l is anomalous. Let us consider 
in the euclidean space the generating function exp[— Z] = J T>ipT>ip exp[— S] ^8|. The most 
general phase transformation onto the Dirac field ip with mass m 

e i(a+&75)0(*) ^ 



induces an additional term in the action given by 

SS = - f d A x [^m (e 2ib ^ e ^ - l) +^(a + bj^d^x) 



(144) 



(145) 



The rotation ( 144 ) gives rise to a nontrivial jacobian due to the noninvariance of the measure 
Dip T>ip. This may be expressed as an additional contribution to the action 



SSi = i / d 4 x6(x) 



( Q ~ b ) Tr F^ V F^ 
8tt 2 » v 



( a + 6 ) TV wW^FW 
87T 2 r *v» 



(146) 



where F^ L ^ U (F( R ^ pv ) is the field strength which couples to the left-handed (right-handed) 
current of the field ip, while F^ v = \e iiVap F ap . Notice that we have absorbed the gauge 
coupling into the definitions of F^ L ' R ^ U and the traces are over the group indices. 
If the rotation ( [144j) corresponds to the baryon number rotation, then 



1 

— and 6 = 0. 
3 



(147) 



Integrating by parts ( |145|) and requiring that the generating function is invariant under the 
baryonic number transformation, we obtain that the baryonic current = J2q ^Q7^Q is 
anomalous 

A/ ( -glF^F^+gfrf^), (148) 



8 T p 



32vr 2 



where Np is the number of fermionic families, F^ v is the field strength of SU{2)l and 
that of U(1)y with coupling constants g<i and g\, respectively, and we have made use of the 
fact that Tr(T a T 6 ) = \5 ab for the SU(2) L generators and of the values Y = 1/6, 2/3, -1/3 
for Ql, ur and dR, respectively. 

Analogously, if we consider the rotation associated to the lepton number 



we obtain 



1 and 6 = 0, 



(149) 



(150) 



where = X^(^7^ + ^£7 M ^) is the leptonic current. The relation ( [L49|) shows explicitly 
that the current associated to B — L is conserved. In fact - since each quark and leton 
family gives the same contribution to the anomaly and each leptonic flavor is conserved in 
the SM at the classical level - the conserved charges are actually three 



l -B 

3 



(151) 



where Lj {i = e, fj,, t) are the leptonic flavors and J2% Li = L. 
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6.2 Topology of SU(2)l and baryon number violation 

In the previous subsection we have described how the chiral anomaly induces the nonconser- 
vation of the baryonic current. We now wish to understand what is the physical significance 
of the terms in the righ-hand side of the Eq. ( j!48| ). First, we note that these terms may be 
reexpressed as 



. N F 

''32^2 



-gjd^ + gld^ 



(152) 



where 



= 2ef»^(&A a p AZ-±g 2 e abB AZA b p AZy 
2e^ (d v B p B a ) , 



and the variation AB of the total baryon number 



B 



d 3 xJ° B 



(153) 



(154) 



in the time interval At is related to the quantity Ncs an d ^cs> called the Chern-Simons 
numbers, in the same time interval 



AB = N F (AN CS - An CS ), 



(155) 



where 



^cs 
^cs 



16vr 2 
16vr 2 



d 3 x 2e ijk TV 
J d 3 x e^ k diBjBk, 



diAjA k + i-g 2 AiAjA k 



where we have defined Ai = A?cr a /2. 

Now, each U(l)y gauge transformation 



(156) 



Bi ^Bi+-{diU Y )Uy\ 

91 



(157) 



with Uy{x) = e taY<yX ^ leaves recs unchanged. On the contrary, there exist SU(2)l gauge 
transformations 



Ai-tUAiU- 1 + -{8iU)U-\ 

92 



which induce a nonvanishing variation of Ncs 

1 



SN CS 



24vr 2 



d 3 xTr 



(diUp-^djU)!/- 1 ^)!/- 1 



ijk 



(158) 



(159) 



This is due to the topological properties of SU{2). The most generic 2x2 unitary ma- 
trix with determinant equal to unity may be expressed as al + ibiicr 1 , with the condition 
a 2 + |6| 2 = 1. Therefore the topology of SU(2) is the same as S 3 , the surface of the hyper- 
sphere in four dimensions (three-sphere) and the gauge transformations are maps from the 
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euclidean space onto SU{2) ~ S* 3 . To clarify this point further, we recall that classically, 
the ground state must correspond to time-independent field configuration with vanishing 
energy density. We have therefore = 0, which means that the field A is a pure gauge, 
A vac = (i/g)(X7U)U~ 1 and we are working in the gauge Aq = 0. Furthermore, we may 
restrict ourselves to the transformations U that have the same limit in all spatial directions. 
We may take this limit to be the identity in the group, U — > 1 as \x\ — ► oo. Under these 
circumstances, all the configurations A vac may be regarded as describing a ground state. 
We have seen that SU(2) is isomorphic to the three-demensional sphere S 3 . On the other 
hand, the whole three-dimensional space with all points at infinity identified is also topo- 
logically equivalent to S 3 . Therefore the gauge transformation U(x) associated with each 
vacuum is a mapping from S 3 onto S 3 . According to the homotopy theory, such mappings 
fall into equivalence classes. Two mappings x — ► U\(x) and x — > E^O?) belong to the same 
class if there exists a continuous transformation from U±(x) to t^O?)- In the case at hand, 
the classes are labeled by positive or negative integer called the winding number. 
We may consider some standard maps 

U<-°\x) = 
U (1 \x) = 

U {n \x) = 

; (i6o) 

It is easy to check that 5Nqs vanishes for J/(°) and any continous transformation of 

U(x) = U (0 \x)(l + ie a {x)a a ), (161) 

where e a (x) — > when \x\ — > oo. On the other hand, <5iVcs does not vanish if we consider 
the transformation 

5N CS = 1. (162) 

The same result is obtained considering continuous deformations of . It is also possible 
to show that - in general - the transformations and their continuous tranformations 
give rise to 

5N CS (V (n) ) = n. (163) 

Therefore, the gauge transformations of SU(2) may be divided in two categories, those 
which do not change the Chern-Simons number, and those which change the Chern-Simons 
number by n, the winding number. 

Let us now consider the SM in the limit in which the mixing angle is zero, i.e. the 
theory is pure gauge SU{2)l theory coupled to the Higgs field <J>. If we choose the gauge 
Aq = 0, we may go from the classical vacuum defined as 

Q^ c = {4 0) = 0.* (0) = (0,v),N cs = 0} , (164) 



1, 

xq + ix • a 



(so + l*1 2 ) 1/2 > 



42 



to an infinite number of other vacua 

= {4 n) = (i/g2)(VUW)(UWy\$W = U^^ \N CS = n) , (165) 

which are classically degenerate and have different Chern-Simons number. Here we have 
denoted the vacuum expectation value (VEV) of the Higgs field by v. 

If we now go back to eq. (|155| ), we are able to understand the connection among 
the baryonic chiral anomaly, the topological structure of SU{2) and the baryon number 
violation. If the system is able to perform a transition from the to the closest 

one Q^lt l \ the Chern-Simons number is changed by one unity and 

AB = AL = N F . (166) 

Each transition creates 9 left-handed quarks (3 color states for each generation) and 3 
left-handed leptons (one per generation) 

3 3 

3EQl + E4~0. (167) 

i=l i=l 



6.3 The sphaleron 



To quantify the probability of transition between two different vacua, it is important to 
understand the properties of the field configurations which interpolate the two vacua and 



"help" the transition. A fundamental result has been obtained in ref. [67], where it was 
found that there exist static configurations (therefore independent from t) which correspond 
to unstable solutions of the equations of motion. These solutions are called sphalerons 
(which in greek stands for "ready to fall") and correspond to saddle points of the energy 
functional and posses Chern-Simons number equal to 1/2. The situation is schematically 
depicted in Fig. 8 

The sphaleron may be identified by considering the minimum energy path among all 
the paths that, in the configuration space, connect two vacua whose Chern-Simons number 
differs by one unit. Along this path, the sphaleron is the configuration of maximum energy 
and is localized is space, even though - contrary to the case of the soliton - is unstable. 

In the limit of vanishing mixing angle, 9\y — > 0, the sphaleron solution has been found 



explicitly by Klinkhamer and Manton [67] and has the following form 



Aidx 1 



- /(saw) U°° d(U°°) 

92 
in 

-j= h(g 2 vr) U<* 



oo\ — 1 




(168) 



where U°°(xi, x 2 , x%) = U^^xq = 0,xi,X2,x^). The energy functional 
E = J d 3 x -FfiFfi + (Vi$)1(Pi$) + V(&) 



(169) 
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Figure 8: Schematical representation of the energy dependence of the gauge configurations 
as a function of the Chern-Simons number. Sphalerons correspond to the maxima of the 
curve. 

may be reexpressed as 



E 



Atxv 

92 



di 



df 



4 37 +7^[/( 1 -/)l 
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92 

g 2 vr, 



dh 



(170) 



where V{4>) is the potential of the Higgs field. 

The functions / and h are the ones which minimize the energy functional ( |170| ) with 
the boundary conditions 



/(0 



CO 



(171) 



and 



1, £ ► CO. 



(172) 



The sphaleron is therefore the solution which interpolates between Q^ c (for £ —* 0) and 
^lac (for £ — > co). The energy and the typical dimensions of the sphaleron configuration are 
basically the result of the competition between the energy of the gauge configuration and 
the energy of the Higgs field. The latter introduces the weak scale into the problem. From 
the quantitative point of view, the potential energy of the Higgs field is less important and 
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for the sphaleron configuration of dimension £, we have 

1 



A 



92f 
Air 



E(AA ~ — (173) 



while the energy of the Higgs field is 

E(<f>) ~ 4irv 2 £. (174) 
Minimizing the sum E(Ai) + E(<f>) we obtain that the typical dimension of the sphaleron is 

4 P ~ 10~ 16 cm, (175) 

92V 

and 

E sp ~ — ~ 10 TeV. (176) 

'92 

A more accurate result may be found by means of variational methods |37| 

E sp = ^b(±), (177) 

92 \92/ 

where B is a function which depends very weakly on \/ g2'- B(0) ~ 1.52 and B(oo) ~ 2.72. 
Including the mixing angle 9w changes the energy of the sphaleron at most of 0.2%. The 
previous computation of the sphaleron energy was performed at zero temeprature. The 
sphaleron at finite temperature - but still in the broken phase - was computed in ]l6|] 
where it was shown that its energy follows approximately the scaling law 

E sp (T) = E sp ^pt, (178) 

where (cf)(T)) is the VEV of the Higgs field at finite temperature in the broken phase. This 
energy may be writen as 

E sp (T) = 2 -^1b(^-), (179) 
aw \92/ 

where m w (T) = ^g 2 ((j)(T)). 

The Chern-Simons number of the sphaleron may be explicitly computing by plugging 
(|168|) into (|156|) p7j and, as mentioned above, one obtains 



^C P S = o- ( 18 °) 



6.4 Baryon number violating transitions 

The probability of baryon number nonconserving processes at zero temperature has been 
computed by 't Hooft [120] and, as we have already mentioned, is highly suppressed by a 
factor exp(— An/aw) ~ 10" 150 . This factor may be interpreted as the probability of making 
a transition from one classical vacuum to the closest one by tunneling, by going through the 
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barrier of ~ 10 TeV corresponding to the sphaleron. An easy way to evaluate this number 
is to remember that the field configurations that describe the transitions (sphalerons or 
instantons in the case of tunneling) are characterized by Ai ~ 1/(52 v) and therefore their 
contribution to the generating function is 

J VA e' s[A] ~ e ~ 1/aw <C 1. (181) 

On the other side one might think that baryon number violating transitions may be obtained 
in physical situations which involve a large number of fields. The contribution to the 
transition amplitude in processes which involve N fields may be estimated as 

/ 1 \ N 

\A B +l\ 2 ~ e -4r/ow = e -4n/ aw -Nloga w _ (182) 

\a w J 

Therefore, if the number of fields involved is about iV ~ l/a\y, the transition probability 
may become of order unity. The sphaleron may be produced by collective and coherent 
excitations containing N > 1/ctw quanta with wavelength of the order of l sp ~ 1/Mw 
At temperatures T ^> Myy, these modes essentially obey statistical mechanics and the 
transition probability may be computed via classical considerations. Note also that at 
temperatures T <C My/ it is no longer possible to deal with classical considerations because 
the Compton wavelength of the thermal excitations ~ T _1 is much larger than the size of 
the sphaleron. 



6.4.1 Baryon number violation below the electroweak phase transition 



After (or during) the electroweak phase transition | 106| by which the SM gauge group 



SU(3)c <8> SU(2)l (g> U(1)y breaks down to U(l) em , the calculation of the baryon number 
violation rate can be done by using the semiclassical approximations ]8q| . The vacuum 
expectation value of the Higgs field (<p(T)) is nonvanishing and the sphaleron configuration 
may be explicitly written down. 

We now want to estimate the transition probability between two different vacua havig 
Chern-Simons number which differ by one unity. One may use an useful analogy. Let us 
consider a pendulum of mass m and be 9 the angle which determines the position of the 
pendulum with respect to the position at rest, 9 = 0. It is clear that the transformation 

6 -> 9 + 2vr, (183) 

may be consider a sort of gauge transformation since the position 9 and 9 + 2ir are indis- 
tinguishable. The periodic potential reads 

V{9) = {mgh){\- cos 9), (184) 

wheer h is the pendulum length. The energy of the corresponding "sphaleron" - the saddle 
point solution to the equation of motion - is V(tt) = 2mgh. According to the classical 
theory, for energies smaller than V(ir), only oscillations around 9 = are possible. How- 
ever, quantum theory predicts a nonvanishing probability of tunneling through the barrier 
separating 9 = from 9 = 2ir, i.e. a complete rotation of the pendulum. 
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Since the solution to the Schrodinger equation reads 



m = Ae -i de'V^m (185) 

the density probability V for penetration from 9 = to = 2tt is 



V ~ |^(27r)| 2 = \A\ 2 e- 2 h* de^2mv{e) (lg6) 

and the quantum tunneling is exponentially suppressed. Imagine now to raise up the tem- 
perature of the system, so that the pendulum coupled to the thermal bath becomes excited 
at higher and higher energies. As the temperature becomes of the order of V(n), it becomes 
possible for the pendulum to reach the position 8 = tt and to roll down to 6 = 2ir. The 
transition rate is therefore 

T(T) oc e - y W/ T , (187) 

and becomes unsuppressed as long as T> V(ir). 

More formally, one has to remember that one of the fundamental objects in statistical 
thermodynamics is the partition function 

Z = Tre~ p6 , (188) 

where H is the Hamiltonian operator and [3 = T _1 . For a scalar field, one may introduce 
the field eingenstates \(ft(x),t) of the Heisenberg picture field operator (ft(x,t) 

4>(x,t)\(ft(x),t)=(ft(x)\(ft(x),t). (189) 

Then the partition function may be written as "summation" over the eigenstates 

Z = ^{(f){x),t = Q\e~ pA \4>{x),t = 0). (190) 

<j>{x) 

We can now make the analogy with the zero temperature case where in the language of 
field theory 

(0"(x),t'V'(5 ? ),t') = (<P"(x),t = 0\e-^ t "- t ^\<j ) >(x),t = 0) 



oc 



J V(pVirexp i dt J d 3 x (ir<j) - H(ir, 0)) 



(191) 



where the path integral is over all the conjugate momenta of 0, n and over all the functions 
satisfying the boundary conditions (ft = (ft"(x) at t" and (ft = (ft'(x) at t'. If, heuristically we 
introduce a variable 

r = it (192) 

and take the limit of integration 

t' = and t" = -i(3, (193) 
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we obtain 

oc J VfiV-K exp J J d 3 x ^in-^ — TC(ir, < 

where now the new boundary conditions are given by 

<f>(P,x) = <p"{x) and 0(0, x) = (j>'(x). 



(194) 



(195) 



By integrating out the conjugate momenta and identifying the boundary conditions, we 
may compute the partition function 







/ V(j) exp 









d 3 x C (cf), d(p) 



where 



and 



1 



i 



(196) 

(197) 
(198) 



The theory at finite temperature may be therefore interpreted as a theory in 3+1 dimensions 
in the euclidean space with periodic boundary conditions on the time coordinate and period 
= 1/T. 

In the limit of very high temperatures only the zero mode of the expansion eft = 
[3~ l Yin e ~* WnT( fr( tjJ n , x), where io n = 2im/f3, is important and the action reduces to —S3/T 
where S3 is euclidean three-dimensional action 



d s x 



(199) 



The transition probability per unit time and unit volume at finite temperature between 
two different minima at (f>\ and (f>2 of a given potential V((f>) for a generic scalar field <j) is 
therefore given at finite temperature by J85| 



V 



A(T) e 



-Sa/T 



(200) 



where A(T) is a prefactor which, on dimensional argument, is 0(T 4 ) and the three-dimensional 
action must be computed for the field configuration (bounce solution) which interpolates 
between the two vacua 

d*l 

dr 2 
lim r _ 



dr 



2d(j) 
r dr 


= v'{4>), 


4>(r) 


= 4>i, 




= 0, 


r=0 




m 


= 02, 


r 


= \x 1. 



(201) 
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This configuration is a bubble whose interior is characterized by the value of the scalar field 
02 and the exterior by 4>i- 



Exercise 3 

Estimate the typical size R c of the bounce solution in the limit of thick bubbles. 

After this long disgression, we are ready to estimate the topological transition rate. 
Since the transition which violates the baryon number is is sustained by the sphaleron 
configuration, one gets S3 = E sp (T). The prefactor was computed in pll] as 



r sp ~ 



2.8 x 10 5 T 4 



fa w \ 4 



B 



-E sp (T)/T 



(202) 



where B has been defined in (170) and k is the functional determinant associated to the 
fluctuations about the sphaleron. It has been estimated to be in the range 10~ 4 < k < 10 _1 



6.4.2 Baryon number violation above the electroweak phase transition 

At temperatures above the electroweak phase transition, the vacuum expectation value of 
the Higgs field is zero, {4>(T)) = 0, the Higgs field decouples and the sphaleron configuration 
ceases to exist. The relevant configuration, at this point, are of the form 



Aidx 1 



0. 
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fU°° d(u°°)- 1 . 



(203) 



Let us estimate the rate r sp on dimensional grounds. As we mentioned, at high temperature 
T the Higgs field decouples from the dynamics and it suffices to consider a pure SU (2) gauge 
theory. Topological transitions take place through the creation of non-perturbative, nearly 
static, magnetic field configurations that generate a change in the Chern-Simons number 
A Acs with a corresponding baryon number generation AB = NjANqs- 

If the field configuration responsible for the transition has a typical scale £, a change 
ANqs — 1 requires 



AA CS 



~ g\ £ 3 dAiA,, 



- g¥ 3 4^ 



This means that the typical energy of the configuration is 



Ai ~ — . 

92i 



(204) 



(205) 



To evade the Boltzmann suppression factor this energy should not be larger than the tem- 
perature T, which requires 



£ > 



1 

gW 



(206) 



49 



Such a length scale corresponds to the one of the dynamically generated magnetic mass of 
order g\T which behaves as a cut off for the maximum coherence length of the system. The 
rate of one unsuppressed transition per volume £ 3 and time t ~ £ is therefore 

T sp ~±- t ~( aw T) 4 . (207) 

This simple scaling argument has been recently criticized in refs. [||, || where it has been 
argued that damping effects in the plasma suppress the rate by an extra power of a\y to 
give r sp ~ a w T 4 . Indeed, since the transition rate involves physics at soft energies g\T 
that are small compared to the typical hard energies ~ T of the thermal excitations in the 
plasma, the simplest way of analyzing the problem is to consider an effective theory for 
the soft modes, where the hard modes have been integrated out and to keep the dominant 



contributions, the so-called hard thermal loops [14|. It is the resulting typical frequency lo c 
of a gauge field configuration immersed in the plasma and spatial with extent (g 2 T)~ 1 that 
determines the change of baryon number per unit time and unit volume. This frequency uo c 
has been estimated to be ~ g%T when taking into account the damping effects of the hard 
modes || |J. This gives r sp ~ j§- ~ a w T 4 . The effective dynamics of soft nonabelian gauge 



fields at finite temperature has been recently addressed also in [13|, where it was found that 
r sp ~ d w T A ln(l/ aw)- Lattice simulations with hard-thermal loops included have been 
performed |)^] and seem to indicate the r sp ~ SOa^T 4 , which is not far from a w T 4 . 
In order to see whether these predictions are reliable, one should write down an effective 
classical hamiltonian for the soft modes of the gauge configurations after having integrated 
out also the soft loops between magnetic fields. These soft loops result to be crucial since 
they not only renormalize the effective aw coupling to non perturbative values, but also 
falsify the naive dimensional arguments about the typical time scale of the sphaleron-like 
fluctuations. From now on, we will parametrize the sphaleron rate as 

T sv = K{a w T)\ (208) 

6.5 The wash-out of B + L 

Let us suppose - for sake of simplicity - that all the charges which are conserved by the 
interactions of the particles in the plasma {Q, Li, B — L, £i = B/3 — Li, ■ ■ •) are zero. If we 
introduce a chemical potential for the charge B + L, hb+l, the free energy density of the 
system (femions) is given by 

F = T I (0 N i 1 + e-^-^' T ) + (jmb+l - -hb+l)] ■ (209) 
The charge density of B + L may be expressed in terms of the chemical potential by 

n B +L ~ Vb+lT 2 (210) 
and - therefore - we may relate the free energy with ub+l 

F ~ ^ B+L T 2 + 0(T 4 ) ~ ^ + 0(T 4 ). (211) 
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The free energy increases quadratically with the fermion number density and the transitions 
which increase ng+L are energetically disfavoured with respect to the ones that decrease 
the fermion number. If these transitions are active for a long enough period of time, the 
system relaxes to the state of minimum energy, i.e. tlb+l = 0: any initial asymmetry in 
B + L relaxes to zero. 

To address this issue more quantitatively, one has to consider the ratio between the 
transitions with 5Nqs = +1 and the ones with 5Nqs = — 1 

^ = e" A ^, (212) 

wheer A/ is the free energy difference between the two vacua. If we define r sp to be the 
average between T+ and T_, we may compute the rate at which the baryon number is 



washed out |lCj 



dn B +L r -r 13 at r s P , 01 „^ 

dt = ~ — — Nf ~7ps n B+L- (213) 

Equation fl213|) is crucial to discuss the fate of the baryon asymmetry generated at the GUT 
scale and is called Master equation. 

Let us now consider temperatures much above the electroweak phase transition, T ^> 
Mw. Baryon number violation processes are active at very high temperatures if the rate 



207) is smaller than the expansion of the Universe 



^>H^T< a tv%~ 10 12 GeV. (214) 

If so, any preexisting asymmetry in B + L is erased exponentially with a typical time scale 
r ~ 2iV>T 3 /13r sp . 

Let us now consider temperatures T ~ Mw when the electroweak phase transition is 
taking place and the Higgs VEV (<j)(T)) is not zero. Baryon number violation processes 
are out-of-equilibrium if, again, the rate ( [20 2D is smaller than the expansion rate of the 
Universe. This translates into the bound on E sp (T) |l0| 

^P>45, (215) 

wheer we have indicated by T c the critical temperature at which the electroweak phase 
transition is taking place. Using the relation ( |179[ ) this bound may be translated into a 
bound on (cp(T c )) 



£1- ( 216 ) 



Any generation of the baryon asymmetry at the electroweak phase transition requires - 
therefore - a strong enough phase transition, that is able to produce a VEV for the Higgs 
field larger than the critical temperature. We will come back to this point later on. 

6.5.1 A crucial point 

In all the considerations leading to Eq. ( |213| ) we have been assuming that all the charges 
which are conserved by the interactions of the particles in the plasma (Q, B — L, Li, 
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ti = B — Li/3, ■ ■ •) are vanishing. Suppose now that these charges - let us denote them 
generically by Qi - are not zero. Define by (B + L)eq the value of the number density 
associated to the B + L charge when the sphaleron transitions are in equilibrium in the 
plasma (ideally, when the sphaleron rate r sp — * oo). In such a case, it is possible to show 
(see Exercise 4) that (B + L)eq is not vanishing in the plasma (6^] 



(B + L) EQ = J2 (kQi, (217) 

where the numerical coefficients Cj depend upon which interactions are in equilibrium in the 
plasma and the particle content of the theory. 

Eq. (217) tells us that anomalous baryon number violating processes do not wash out 



completely the combination B + L if at least one of the charges which are conserved by the 
interactions of the plasma, e.g. (B — L) is nonvanishing. Eq. ( [21 3D changes accordingly 
(see Exercise 4): 



dn B+ L r sp dF r sp 

DC DC — 

dt T 8(B + L) T 



EQ 

HB+L) ~ n (B+L) 



(218) 



The interpretation of this equation is straightforward. It is a Boltzmann equation in the 
sense that, in the limit r sp —* oo, the solution is B + L = [B + L)eq- If the conserved 
charges are zero, then any B + L is washed-out, see Eq. (|217j ). However, if (B + L)eq 
is not zero, then sphalerons transitions will act on the system until the (B + L) charge 
has been reduced to its equilibrium value. The latter is not necessarily zero if some other 
conserved charge, like B — L, is not zero. In other words, sphaleron transitions push the 
system towards the state of minimum free energy, which is characterized by a nonvanishing 
B + L if other conserved charges are non zero. This is a crucial point to keep in mind when 
we will talk about electroweak baryogenesis. 



Exercise 4 



a) Consider the two Higgs doublet model in the broken phase. The Higgs doublets are 
defined as Hi = (H®, H~) T and H2 = (H + , iJ^ ) T an d couple to the down-type and up-type 
quarks, respectively. By considering all the processes in thermal equilibrium (but the ones 
mediated by light quark Yukawa interactions, Cabibbo suppressed gauge interactions and 
sphalerons transitions), identify the charges which are conserved by the interactions. Hint: 
one of them is B + L; b) assuming that also the sphaleron transitions are in equilibrium, 
compute the relation between the corresponding equilibrium value of the B + L charge (call 
it (B + L)eq), and the other conserved charges, c) Compute the free energy of the system 
and show that it scales like [(B + L) — (B + L)eq] 2 . 

6.6 Baryon number violation within the SM and GUT baryogenesis 

At this point, we are ready to discuss the implications of the baryon number violation in 
the early Universe for the baryogenesis scenarios discussed so far. The basic lesson we 
have learned in the previous subsections is that any asymmetry B + L is rapidly erased by 
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sphaleron transitions as soon as the temperatures drops down ~ 10 12 GeV. Now, we can 
always write the baryon number B as 

B + L B-L , , 

B = — + —- ^ 

This equation seems trivial, but is dense of physical significance! Sphaleron transitions only 
erase the combination B + L, but leave untouched the orthogonal combination B — L. This 
means that the only chance for a GUT baryogenesis scenario to work is to produce at high 
scale an asymmetry in B — L. In section 4 - however - we have learned that there is no 
possibility of generating such an asymmetry in the framework of SU(5). This is because 
the fermionic content of the theory is the one of the SM and there is no violation of B — L. 
Sphaleron transitions are therefore the killers of any GUT baryogenesis model based on 
the supersymmetric version of SU(5) with R parity conserved (The non-supersymmetric 
version is already ruled out by experiments on the proton decay lifetime). This is a striking 
result. 

6.6.1 Baryogenesis via leptogenesis 

The fact that the combination B — L is left unchanged by sphaleron transitions opens up 
the possibility of generating the baryon asymmetry from a lepton asymmetry. This was 
suggested by Fukugita and Yanagida p7[| . The basic idea is that, if an asymmetry in the 
lepton number is produced, sphaleron transition will reprocess it and convert (a fraction of) 
it into baryon number. This is because B + L must be vanishing all the times and therefore 
the final baryon asymmetry results to be B ~ —L. The primordial lepton asymmetry 
is generated by the out-of-equilibrium decay of heavy right-handed Major ana neutrinos 
N L . Once the lepton number is produced, the processes in thermal equilibrium distribute 
the charges in such a way that in the high temperature phase of the standard model the 
asymmetries of baryon number B and of B — L are proportional in thermal equilibrium |62|] 
(see Eq. (Q 

Where Njj is the number of Higgs doublets. As we have already stressed, in the standard 
model, as well as its unified extension based on the group SU(5), B — L is conserved. 
Hence, no asymmetry in B — L can be generated, and B vanishes. However, a nonvanishing 
B — L asymmetry may be naturally obtained adding right-handed Majorana neutrinos to 
the standard model. This extension of the standard model can be embedded into GUTs 
with gauge groups containing SO(W). Heavy right-handed Majorana neutrinos can also 



explain the smallness of the light neutrino masses via the see-saw mechanism [50]. 

The basic piece of the Lagrangian that we need to understand leptogenesis is the coupling 
between the right-handed neutrino, the Higgs doublet <I> and the lepton doublet II 

£ = I£$h„N£ + -N£M N£ + h.c. (221) 

The vacuum expectation value of the Higgs field (#) generates Dirac masses for neu- 
trinos mo = h v {$), which are assumed to be much smaller than the Majorana masses 
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<♦> 

(b) 



Figure 9: (a) One-loop interference giving rise to the lepton asymmetry; (b) Diagram giving 
rise to the 5-dimensional operator of Eq. (216). 



M . When the Major ana right-handed neutrinos decay into leptons and Higgs scalars, they 
violate the lepton number ( right-handed neutrino fermionic lines do not have any preferred 
arrow) 



$ + 
$+1 



(222) 



The interference between the tree- level and the one-loop amplitudes, see Fig. 9(a), yields a 
CP asymmetry equal to 



where 



87ru 2 (jn^mnj 



m} D mD 



i.i 



1 \M?. 



1 - (1 + x) In 



1 + x 



(223) 



(224) 



and the index i is summed over all the three species of right-handed neutrino. The final 
baryon asymmetry has been computed by several authors |9l], 81, 3C, 17] and it has been 
shown to be of the order of 

B ~ (0.6 - 1) x 1(T 10 . (225) 
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However, one has to avoid a large lepton number violation at intermediate temperatures 
which may potentially dissipate away the baryon number in combination with the sphaleron 
transitions. Indeed, the diagram of Fig. 9(b), induced by the exchange of a heavy right- 
handed neutrino, gives rise to a AL = 2 interaction of the form 

v £ L £ L $$ + h.c, (226) 



where m u is the mass of the light left-handed neutrino. The rate of lepton number violation 
induced by thi interaction is therefore Tl ~ (m 2 /(<l?) 4 )T 3 . The requirement of harmless 
letpon number violation, <H imposes an interesting bound on the neutrino mass 

'"^ev( I ^)~ 1/2 , (227) 

where Tx = Min {Tb_l, 10 12 GeV} and Tb~l is the temperature at which the B—L number 
production takes place and ~ 10 12 GeV is the temperature at which sphaleron transitions 
enter in equilibrium. One can also reverse the argument and study leptogenesis assuming 
a similar pattern of mixings and masses for leptons and quarks, as suggested by 50(10) 
unification [17|. This implies that B — L is broken at the unification scale ~ 10 16 GeV, 



if m u ~ 3 x 10 eV as preferred by the MSW explanation of the solar neutrino deficit 



[ 1251 , |7|. 



7 Electroweak baryogenesis 

So far, we have been assuming that the departure from thermal equilibrium, necessary 
to generate any baryon asymmetry, is attained by late decays of heavy particles. In this 
section, we will focus on a different mechanism, namely the departure from equilibrium 
during first order phase transitions. 

A first order phase transition is defined to occur if some thermodynamic quantities 
change discontinuously. This happens because there exist two separate thermodynamic 
states that are in thermal equilibrium at the time of the phase transition. The thermody- 
namic quantity that undergoes such a discontinuous change is generically called the order 
parameter cf). Whether a phase transition is of the first order or not depends upon the 
parameters of the theory and it may happen that, changing those parameters, the order 
parameter becomes continuous at the time of the transition. In this case, the latter is said 
to be of the second order at the point at which the transition becomes continuous and a con- 
tinuous crossover at the other points for which all physical quantities undergo no changes. 
In general, we are interested in systems for which the high temperature ground state of the 



theory is at <j) = and the low temperature phase is at 4> / |7^, 106 1- 

For a first order phase transition, the extremum at <p = becomes separated from a 
second local minimum of the potential by an energy barrier. At the critical temperature 
T c both phases are equally favoured energetically and at later times the minimum at <j) ^ 
becomes the global minimum of the theory. The phase transition proceeds by nucleation 
of bubbles. Initially, the bubbles are not large enough for their volume energy to overcome 
the competing surface tension, they shrink and disappear. However, at the nucleation 
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Figure 10: Schematic picture of the propagating bubble separating the broken from the 
unbroken phase during the electroweak phase transition. 

temperature, critical bubbles form, i.e. bubbles which are just large enough to be nucleated 
and to grow. As the bubble walls separating the broken from the unbroken phase pass each 
point in space, the order parameter changes rapidly, leading to a significant departure from 
thermal equilibrium. 

The critical bubbles of the broken (Higgs) phase have a typical profile 



where r is th spatial coordinate, is the bubble wall width and {<p{T c )) is the VEV of the 
Higgs field inside the bubble. 

Bubbles expand with velocity until the fill the Universe; local departure from thermal 
equilibrium takes place in the vicinity of the expanding bubble walls, see Fig. 10. 

The fundamental idea of electroweak baryogenesis is to produce asymmetries in some 
local charges which are (approximately) conserved by the interactions inside the bubble 
walls, where local departure from thermal equilibrium is attained. These local charges will 
then diffuse into the unbroken phase where baryon number violation is active thanks to 
the unsuppressed sphaleron transitions. The latter convert the asymmetries into baryon 
asymmetry, because the state of minimum free energy is attained for nonvanishing baryon 
number, see eq. ( |218| ). Finally, the baryon number flows into the broken phase where it 
remains as a remnant of the electroweak phase transition if the sphaleron transitions are 
suppressed in the broken phase. The recipe for electroweak baryogenesis is therefore the 
following: 




(228) 
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- Look for those charges which are approximately conserved in the symmetric phase, 
so that they can efficiently diffuse in front of the bubble where baryon number violation is 
fast, and non-orthogonal to baryon number, so that the generation of a non-zero baryon 
charge is energetically favoured. 

- Compute the CP violating currents of the plasma locally induced by the passage of 
the bubble wall. 

- Write and solve a set of coupled differential diffusion equations for the local particle 
densities, including the CP violating source terms derived from the computation of the 
current at the previous step and the particle number changing reactions. The solution to 
these equations gives a net baryon number which is produced in the symmetric phase and 
then transmitted into the interior of the bubbles of the broken phase, where it is not wiped 
out if the first transition is strong enough. 



7.1 Electoweak baryogenesis in the SM 



Since C and CP are known to be violated by the electroweak interactions, it is possible 
- in principle -to satisfy all Sakharov's conditions within the SM if the electroweak phase 
transition leading to the breaking of SU(2)l <S> U(1)y is of the first order [62]. There are 
very good reviews on electroweak baryogenesis and the reader is referred to them for more 
details | 



113, 122, 38]. 



The asymmetry flowing inside the bubbles of the broken phase will survive if sphaleron 
transitions are frozen out and baryon number violation is inefficient. As we have learned in 
the previous section, baryon number violation is out-of-equilibrium inside the bubble wall 
only of > 1, i.e. if the electoweak phase transition is strong first order. Let us now 

understand as this condition translates into a upper bound on the Higgs mass m^. 

In general, given an order parameter 4> an d a set of particles i with masses rnj(^) in 
the <j> background, plasma masses vrj(T) and degrees of freedom m, the effective one- loop 
improved potential at finite temperature is given by [fU06 



AV 



bos / 



■>,t) = Y. 



I 24 12vr 



m^) + IT(T) 



3/2 mt(6) , mf{6) ] . 

log -p-T^r \ (229) 



64vr 2 



AbT 2 j 



if the particles are bosons and 



AV 



for / 
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(230) 



if they are fermions. Here Ab = 16 Ap = 16-k 2 exp(3/2 — 2je), 7e — 0.5722. 

One can therefore write the total one-loop effective potential of the SM Higgs field at 



finite temperature as as [ 106 ] 



V(cf), T) = D(T 2 - T$)4) 2 - ETct> 3 + ^-V, 



(231) 



where 



D 



2M' 2 V + M| + 2m 2 



8v 2 
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Id ' 



A(T) = A - ^ ( 2 M^ Log£f, + Mi Log-M,- 4^ Log^) , (232) 

where mj is the mass of the top-quark. 

It is now easy to see that, when the minimum (j> = becomes metastable, i.e at the 
temperature T c when V(0, T c ) = V((p(T c ),T c ), one has 

*P>=lfU^. (233) 

where we have used the fact that m| = 2Ai> 2 . The condition fl216| ) is therefore satisfied only 
if 

[ae 

m h <\ 42 GeV. (234) 

y 1.3 

On the other hand, the current lower bound on comes from combining the results 



of DELPHI, L3 and OPAL experiments and is > 89.3 GeV [11]. A simple one-loop 
computation shows, therefore, that the electroweak phase transition is too weakly first 
order to assure the preservation of the generated baryon asymmetry at the electroweak 
phase transition in the SM. More complete perturbative and non-perturbative analyses 
|113| 1 have shown that the electroweak phase transition is first order if the mass of the Higgs 
nth is smaller than about 80 GeV and for larger masses becomes a smooth crossover. Let 
us now briefly analyzed the issue of CP violation within the SM. Because of CP violation 
in the kaon system, it is of great interest to see whether enough CP violation is present in 
the SM to generate the baryon asymmetry at the observed level. 

A very rough (and optimistic) estimate of the amount of CP violation necessary to 
generate B ~ 10 -10 can be obtained as follows. Since the baryon number violation rate in 

factor due to CP violation, we get 



the symmetric phase is proportional to aA, ~ 10 6 , if we indicate by 5c p the suppression 



rv 4 T 3 

B ~ 5 C p ^ 10~ 8 5 C p. (235) 

s 

Even neglecting all the suppression factors coming from the dynamics of the electroweak 
phase transition, we discover that 

&CP >10~ 3 . (236) 

A naive estimate suggests that, since CP violation vanishes in the SM if any two quarks of 
the same charge have the same mass, the measure of CP violation should be the Jarlskog 
invariant 

^f- = (m?-mZ)(m 2 c -mI)(m 2 u -m?) 

(M 2 - M 2 ) (M 2 - M 2 ) (M 2 - M 2 ) , (237) 
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where J is twice the area of the unitarity triangle. The quantity Aqp has dimension twelve. 
In the limit of high temperature, T much larger than the quark masses M q , the only mass 
scale in the problem is the temperature itself. Therefore, the dimensionless quantity 5c p is 

SCP * * 10- 20 , (238) 

far too small for the SM to explain the observed baryon asymmetry. 

This admittedly too naive reasoning has been questioned by Farrar and Shaposhnikov 



[45] who have pointed out that for quarks having momentum p ~ T 3> M q , the above 
estimate is certainly correct since light quarks are effectively degenerate in mass and the 
GIM suppression is operative; on the other side, this is no longer true when quarks have a 
momentum p ~ M q . Since the mass jump through the bubble wall is just M q , quarks coming 
from the symmetric phase and with momentum p < M q are reflected off from the wall, while 
the ones with momentum p > M q are partially reflected and partially trasmitted. In the 
reflection processes quarks and antiquarks acquire different probabilities of penetrating the 
bubble wall. In such a way, it might be possible to produce a net baryon number flux from 
outside to inside the bubble wall. For instance, considering momenta between and M s , 
then all the strange quarks might be reflected off, while down quarks have a nonvanishing 
probability of being transmitted. However, this effect is largely suppressed by the fact that 
fermions, when they propagate in the plasma, acquire a damping rate 7 ~ 0.1 T 3> M s and 
the quark energy and momenta cannot be defined exactly, but have a spread of the order of 
7 S> {M a — Mfj). In other words, the lifetime of the quantum packet is much shorter than 
the typical reflection time from the bubble wall (~ 1/M S ): CP violation, which is based 



on coherence and needs at least a time ~ 1/M S to be built up, cannot be efficient [48, 59]. 
Therefore, the common wisdom is that electroweak baryogenesis is not possible within the 
SM. 

7.2 Electoweak baryogenesis in the MSSM 

The most promising and well-motivated framework for electroweak baryogenesis beyond the 
SM seems to be supersymmetry (SUSY) |56| , |[] . Let us remind the reader only a few notions 
about the MSSM that will turn out to be useful in the following. 
Let us consider the MSSM superpotential 

W = fiHiH 2 + h u H 2 Qu c + h d HiQd c + h e H x Le c , (239) 

where we have omitted the generation indices. The Higg sector contains the two Higgs 
doublets 

H, = j and H 2 =^y (240) 

The lepton Yukawa matrix h e can be always taken real and diagonal while h u and h d contain 
the KM phase. 

What is relevant for baryogenesis is to identify possible new sources of CP violation. 
They emerge from the operators which break softly supersymmetry 
i) Trilinear couplings: 
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T u H 2 Qu c + ^HiQd + V e H x L^ + h.c, 



(241) 
(242) 



where we have defined 

T (u,d,e) = ^ A (u,d,e) . h (u,d,e)_ 

Generally, in supergravity models the matrices are assumed to be proportional to 

the identiy matrix 

A^' d ' e \M GVT ) = A-l, (243) 



where the A parameter can be complex. 
ii) bilinear couplings: 



fiBHiH 2 + h.c. 



Hi) Majorana gaugino masses: 
1 



(M1A1A1 + M 2 A 2 A 2 + M3A3A3) + h.c. 
At the GUT scale it is usually assumed that 

Mi = M 2 = M 3 = M. 

iv) Scalar soft masses: 

ml b z* a z + h.c. 



(244) 
(245) 
(246) 

(247) 



The new contributions to explicit violation of CP are given in the phases of the complex 
parameters A, B, Mi (i = 1, 2, 3) and by the parameter [i in the superpotential ( |239| ) . Two 
phases may be removed by redifining the phase of the superfield H 2 in such a way that the 
phase of \i is opposite to that of B. The product fiB in ( 244[ ) is therefore real. It is also 
possible to remove the phase of the gaugino mass M by an R symmetry transformation. 
The latter leaves all the other supersymmeric couplings invariant and only modifies the 
trilinear ones, which get multiplied by exp(— 4>m) where 4>m is the phase of M. 
The phases which are left are therefore 



arg(AM) and 



-arg(.B). 



(248) 



The two new phases 4>a and (j)^ will be crucial for the generation of the baryon asymmetry. 

Electroweak baryogenesis in the framework of the Minimal Supersymmetric Standard 
Model (MSSM) has attracted much attention in the past years, with particular emphasis on 



the strength of the phase transition [51, |100| , 42, 15] and the mechanism of baryon number 
generation 1 102, 




Recent analytical [fLSj,~l3*3|i 



9[ and lattice computations J77|, |27, 



78, [7{|] have revealed that the phase transition can be sufficiently strongly first order if the 
ratio of the vacuum expectation values of the two neutral Higgses tan (3 is smaller than ~ 4. 
Moreover, taking into account all the experimental bounds as well as those coming from the 
requirement of avoiding dangerous color breaking minima, the lightest Higgs boson should 
be lighter than about 105 GeV, while the right-handed stop mass might be close to the 
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present experimental bound and should be smaller than, or of the order of, the top quark 
mass [p0| . 

Moreover, as we have seen, the MSSM contains additional sources of CP-violation be- 
sides the CKM matrix phase. These new phases are essential for the generation of the 
baryon number since large CP violating sources may be locally induced by the passage 
of the bubble wall separating the broken from the unbroken phase during the electroweak 
phase transition. Baryogenesis is fuelled when transport properties allow the CP violating 
charges to efficiently diffuse in front of the advancing bubble wall where anomalous elec- 
troweak baryon violating processes are not suppressed. The new phases appear in the soft 
supersymmetry breaking parameters associated to the stop mixing angle and to the gaugino 
and neutralino mass matrices; large values of the stop mixing angle are, however, strongly 
restricted in order to preserve a sufficiently strong first order electroweak phase transition. 
Therefore, an acceptable baryon asymmetry from the stop sector may only be generated 
through a delicate balance between the values of the different soft supersymmetry breaking 
parameters contributing to the stop mixing parameter, and their associated CP violating 



phases [19|. As a result, the contribution to the final baryon asymmetry from the stop sector 



turns out to be negligible. On the other hand, charginos and neutralinos may be responsible 



for the observed baryon asymmetry if the phase of the parameter // is large enough [19, 28]. 
Yet, this is true within the MSSM. If the strength of the electroweak phase transition is 
enhanced by the presence of some new degrees of freedom beyond the ones contained in 
the MSSM, e.g. some extra standard model gauge singlets, light stops (predominantly the 
right-handed ones) and charginos/neutralinos are expected to give quantitatively the same 
contribution to the final baryon asymmetry. 

7.2.1 The electroweak phase transition in the MSSM 

As discussed above, a strongly first order electroweak phase transition can be achieved in the 
presence of a top squark lighter than the top quark |20]]. In order to naturally suppress its 
contribution to the parameter Ap and hence preserve a good agreement with the precision 
measurements at LEP, it should be mainly right handed. This can be achieved if the left 
handed stop soft supersymmetry breaking mass itiq is much larger than Mz- 
The stop mass matrix is given by 

( S "k ) ■ (249) 



where 



M 2 LL * m 2 Q + h 2 t 
M 2 RR ^ ml + h 2 



Hi 



M 2 LR = h t (A t H%-ii*Hl). (250) 
For moderate mixing, the lightest stop mass is then approximately given by 

A 2 " 

m 2 7 ~mtj + mi{(t))\l | (251) 
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where At = At — [i*/ tan (5 is the particular combination appearing in the off-diagonal terms 
of the left-right stop squared mass matrix and mfj is the soft supersymmetry breaking 
squared mass parameter of the right handed stop. Notice that the Higgs sector contains 
two neutral CP even states, and H®. However, in the limit in which itia 3> T c , where 
rriA is the mass of the pseudoscalar particle of the Higg sector, only one neutral Higgs 
survives 

<j> = cos PH$ + sin PH$, (252) 

where tan/3 = {H$)/(Hf}, and the low -energy potential reduces to the one-dimensional 
SM-like potential V (</>). 

The preservation of the baryon number asymmetry requires the order parameter (<f>(T c )) /T c 
to be larger than one. The latter is bounded from above 



<*(K» < (mL) + w(i-i/r„ t (253) 



3/2 



T c \ T c J SM irvm 2 h 

where mt = rrit(mt) is the on-shell running top quark mass in the MS scheme. The first 
term on the right hand side of expression ( 253| ) is the Standard Model contribution 
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( (<KT C )) 
V T c Jsm \m h [GeV] 



(254) 



and the second term is the contribution that would be obtained if the right handed stop 
plasma mass vanished at the critical temperature (see Eq. ( |255| )). Remember that in the 
expression for the one-loop effective potential fl231| ), the parameter E gets contributions 
from boson fields. So, the difference between the SM and the MSSM is that light stops may 
give a large contributions to the effective potential in the MSSM. 

In order to overcome the Standard Model constraints, the stop contribution must be 
therefore large. The stop contribution strongly depends on the value of mfj, which must be 
small in magnitude, and negative, in order to induce a sufficiently strong first order phase 
transition. Indeed, large stop contributions are always associated with small values of the 
right handed stop plasma mass 

mf = -ml + Ii R {T), (255) 

where rhfj = —mfj, Hr(T) ~ %|T 2 /9 + /tf/6[2 — Af/mfyT 2 is the finite temperature self- 
energy contribution to the right-handed squarks. Moreover, the trilinear mass term, At, 
must be A\ <C rriq in order to avoid the suppression of the stop contribution to (<p(T c )) /T c . 

Although large values of rhu, of order of the critical temperature, are useful to get a 
strongly first order phase transition, they may also induce charge and color breaking minima. 
Indeed, if the effective plasma mass at the critical temperature vanished, the universe would 
be driven to a charge and color breaking minimum at T > T c . Hence, the upper bound 
on (</>(T c ))/T c , Eq. Q253] ) cannot be reached in realistic scenarios. A conservative bound 
on fhjj may be obtained by demanding that the electroweak symmetry breaking minimum 
should be lower than any color-breaking minima induced by the presence of fhjj at zero 
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temperature, which yields the condition 



mu < 




(256) 



It can be shown that this condition is sufficient to prevent dangerous color breaking minima 
at zero and finite temperature for any value of the mixing parameter At- A more general 
analysis is provided in [pCfl . 

In order to obtain values of (c/>(T c )) /T c larger than one, the Higgs mass must take small 
values, close to the present experimental bound. Numerically, an upper bound, of order 80 
GeV, can be derived. For small mixing, the one-loop Higgs mass has a very simple form 



where m~ ~ uiq + mf, is the heaviest stop squared mass. Hence, tan/3 must take values 
close to one. The larger the left handed stop mass, the closer to one tan/3 must be. This 
implies that the left handed stop effects decouple at the critical temperature and hence, 
different values of uiq mainly affect the baryon asymmetry through the resulting Higgs 
mass. 

Values of the CP-odd Higgs mass tjia 200 GeV are associated with a weaker first 
order phase transition. Fig. 11 shows the behaviour of the order parameter {(p(T c )) /T c in 
the m^-tan (3 plane, for A t = 0, rriQ = 500 GeV and values of fhy close to its upper bound, 



In order to correctly interpret the results of Fig. 11 one should remember that the Higgs 
mass bounds are somewhat weaker for values of rriA 150 GeV. However, even for values 
of rriA of order 80 GeV, in the low tan (5 regime the lower bound on the Higgs mass is of 
order 60 GeV. Hence, it follows from Fig. 11 that, to obtain a sufficiently strong first order 
phase transition the CP-odd Higgs mass txia ^ 150 GeV. When two-loop QCD corrections 
|^3| , |44j associated with stop loops are included, one finds that tua > 120 GeV irih < 85 
GeV. This region will be explored at LEP2 very soon. 



Obtain Eq. (|253; ). 



7.2.2 How to produce the baryon asymmetry in the MSSM 

As we have previously learned, the first step in the computation of the baryon number 
asymmetry is to identify those charges which are approximately conserved in the symmetric 
phase, so that they can efficiently diffuse in front of the bubble where baryon number 
violation is fast, and non-orthogonal to baryon number, so that the generation of a non- 
zero baryon charge is energetically favoured according to the Master equation (|218j ). 

Charges with these characteristics in the MSSM are the axial stop charge and the Hig- 
gsino charge, which may be produced from the interactions of squarks and charginos and/or 




(257) 



Eq. (|H). 



Exercise 5 
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100 150 200 250 300 



m A (GeV) 

Figure 11: Contour plots of constant values of (4>(T C ))/T C (solid lines) and in GeV 
(dashed lines) in the plane (mA, tan (3). We have fixed mt = 175 GeV and the values of 
syper symmetric parameters: rriQ = 500 GeV, mjj = m™ 1 fixed by the charge and color 
breaking constraint, and At = ji/ tan/3. 
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neutralinos with the bubble wall, provided a source of CP-violation is present in these sec- 
tors. This is exactly the case, since both the parameters At and \i may carry a physical 
phase 0. The idea is that, if nonvanishing CP violating sources for the right-handed stop 
and higgsino numbers are induced in the bubble wall, the scattering among particles as 
well as diffusion will generate an asymmetry in the left-handed fermion asymmetries in 
the unbroken phase. The asymmetry - in turn - will fuel baryogenesis because sphaleron 
transitions will push the system towards the state of minimum free energy, which is the 
one with nonvanishing baryon asymmetry. In the next subsection, we will give some indi- 
cations of how to compute the CP-violating sources. Let us now investigate the dynamics 
of electroweak baryogenesis a little bit further. 

One has to start with a set of coupled differential equations describing the effects of 
diffusion, particle number changing reactions and CP-violating source terms. Major sim- 
plifications of the diffusion equations take place when neglecting all the couplings except 
for gauge interactions and the top Yukawa coupling. Neglecting the weak sphalerons (in 
the first step) allows to forget about leptons in the diffusion equations and will turn out to 
be a good approximation when computing Higgs and quark densities. 

If the system is near thermal equilibrium and particles interact weakly, the particle 



number densities rij may be expressed as (see Eq. (11)) rii = /c,^jT 2 /6 where \ii is the 
local chemical potential, and fcj are statistical factors of the order of 2 (1) for light bosons 
(fermions) in thermal equilibrium, and Boltzmann suppressed for particles heavier than T. 

What really determines which are the interactions in equilibrium is the typical time 
scale for the passage of the bubble wall through a given point, r w ~ L^/v^. If interactions 
are faster than they are in equilibrium, otherwise, they are not. 

The particle densities we need to include are 

- the left-handed top doublet q = (tz + 

- the right-handed top quark t = tji, 

- the Higgs particle h = (H®, H® , il-f, H%), and the superpartners q, t and h. 
The interactions able to change the particle numbers are 

- the top Yukawa interaction with rate Tt, 

- the top quark mass interaction with rate T m , 

- the Higgs self- interactions in the broken phase with rate I\, 

- the strong sphaleron interactions with rate r ss . The axial vector current of QCD 
J2i ^7^75(7* = — J2i QlI^Ql + Si QrI^Qr where the sum is over the quarks, has a triangle 
anomaly and therefore one may expect axial charge violation due to topological transitions 
analogous to the case of sphaleron transitions 

dQ 5 _ 12-6 

where Q5 is the axial charge, the factor 12 comes from the total number of quark chirality 
states and the factor 6 from the relation between the asymmetry in the quark number 
density and the chemical potential, rtj ~ ^T 2 /6, see Eq. (p"T|). The rate of these processes 



at high temperature is expected to be |92j 



8 /as\ 
3 \a w ) 



r sp , (259) 
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where as is the strong fine structure leading to the characteristic time of order of 

' <t w . (260) 



192 Ka s T 

The effect of QCD sphalerons may be represented by the operator 

nf =1 («i44 4)i. ( 261 ) 

where i is the generation index. Assuming that these processes are in equilibrium |95|, (52J, 
we get the following equation for the chemical potentials 

3 

-Mujj + Z^ — Mdjj) = - ( 262 ) 

i=i 

This equation contains the chemical potential for all the quarks and imposes that the total 
right-handed baryon number is equal to the total left-handed one. In other words, including 
strong QCD sphalerons allow the generation of the right-handed bottom quark as well as 
the generation of the first and second family quarks, 

- the weak anomalous interactions with rate r sp , 

- the gauge interactions. 

We shall assume that the supergauge interactions are in equilibrium, that is 

q q t t h h . 

(263) 



kg k~ kt k t kh 



Under these assumptions the system may be described by the densities Q = q + q, 
T = t + t and H = h + h. CP- violating interactions with the advancing bubble wall 
produce source terms 7^ for Higgsinos and 7^ for right-handed stops, which tend to push 
the system out of equilibrium. Ignoring the curvature of the bubble wall, any quantity 
becomes a function of the coordinate z = Z3 + v u z, the coordinate normal to the wall 
surface, where we assume the bubble wall is moving along the Z3-axis. 

When including the strong sphalerons, right-handed bottom quarks are generated as 
well as the quarks of the first two families. However, since strong sphalerons are the only 
processes which produce the first two generation quarks and all quarks have nearly the same 
diffusion constants, we may constrain the densities algebrically in terms of B = bn + 6r 

Qil = Q2L = -2U R = -2D R = -2S R = -2C R = -2B = 2(Q + T), (264) 

where the last equality comes from imposing that strong sphalerons are in equilibrium. 

Particle transport is treated by including a diffusion term. Taking all the quarks and 
squarks with the same diffusion constant D q and Higgs and Higgsinos with diffusion constant 
Dh, onecan write the folowing set of diffusion equations 

Q = Dq V 2 Q-r t [Q/k Q -H/k H -T/k T ]-r m [Q/k Q -T/k T ] 

- 6T ss [2Q/k Q -T/k T + 9{Q + T)/k B }+ 1I , 

f = D q V 2 T-T t [-Q/k Q + H/k H + T/k T )-T m [-Q/k Q + T/k T ) 

+ 3T ss [2Q/k Q -T/k T + 9{Q + T)/k B }+ 1I , 

H = D h V 2 h-T t [-Q/k Q + H/k H + T/k T ]-T h H/k H + ll , (265) 
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where we have inserted the CP violating sources. 

Assuming that the rates Tt and r ss are fast so that Q/k q — H/kn — T/k? = 0(1/Tf) and 
2Q/k q — T/kx + 9(Q + T)/kb = 0(1/T SS ), one can find the equation governing the Higgs 
density 

v u H' - ~DH" + FH — 7 = 0, (266) 

where the derivatives are now with respect to z, D is the effective diffusion constant, 7 
is an effective source term in the frame of the bubble wall and T is the effective decay 
constant [ |102| ] . An analytical solution to Eq. ( |266| ) satisfying the boundary conditions 
H(±oo) = may be found in the symmetric phase (defined by z < 0) using a z-independent 
effective diffusion constant and a step function for the effective decay rate T = T9(z). A 
more realistic form of T would interpolate smoothly between the symmetric and the broken 



phase values. The values of D and V in (266) of course depend on the particular values of 



supersymmetric parameters. For the considered range one typically finds D ~ 0.8 GeV 1 , 
r ~ 1.7 GeV. 

The tunneling processes from the symmetric phase to the true minimum in the first 
order phase transition of the Higgs field in the MSSM has been recently analyzed in j9t| 
including the leading two-loop effects. It was shown that the Higgs profile along the bubbles 
at the time when the latter are formed has a typical thickness L w ~ (20 — 30) /T. In 
general, however, the value of when the bubbles are moving through the plasma with 
some velocity v w is different from the value at bubble nucleation. Indeed, the motion of 
the bubble wall is determined by two main factors, namely the pressure difference between 
inside and outside the bubble -leading to the expansion- and the friction force, proportional 
to v u , accounting for the collisions of the plasma particles off the wall. The equilibrium 
between these two forces imples a steady state with a final velocity v w . If bubbles are rather 
thick, thermodinamical conditions are established inside the wall and for the latter is no 
longer possible to loose energy by thermal dissipation. Under these conditions the bubble 
wall is accelerated until slightly out-of-equilibrium conditions and the friction forces are 
reestablished. As we shall see, the total amount of the baryon asymmetry is proportional to 
A/3 -the change in the ratio of the Higgs vacuum expectation values j3 = {H®) / {H®) from 
z = to inside the bubble wall. This quantity tends to zero for large values of itia, and 
takes small values, of order 10 -2 for values of the pseudoscalar mass = 150-200 GeV 



The solution of Eq. ( 266 ) for z < is 

H(z) =Ae zv "/~ 5 , (267) 

and for z > is 

* w = ( fi+ -w^r^ (u)e " + ") ei+z 

+ 1 - - f Z duj{u)e~ x - u ] e x - z . (268) 

where 



L>(A_ - A+) Jo 



= ^ — 269 

2D y 1 
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Imposing the continuity of H and H' at the boundaries, we find 



A_\ „ f\+ A 1 



A = B +{ 1 -r + )= B -{T-- 1 ) = Tf>r + L ^y~^- < 2 ™> 

From the form of the above equations one can see that CP violating densities are non zero 
for a time t ~ D/v^ and the assumptions leading to the analytical form of H(z) are valid 
provided that the interaction rates and r ss are larger than vl/D g(j|, 
The equation governing the baryon asymmetry n B is given by [102] 



D q n" B - v u n' B - 9(-z)N f T sp n L = 0, (271) 

where ni is the total number density of left-handed weak doublet fermions and we have 
assumed that the baryon asymmetry gets produced only in the symmetric phase. Expressing 
ra^(z) in terms of the Higgs number density 

9k q k T - 8k b k T - 5k b k q 
nL ~ k H (k b + 9k q + 9k T ) H 

and making use of Eqs. (|267| )- (|27lD , we find that 

2» = - g ( h )^M, (273) 

where g{ki) is a numerical coefficient depending upon the light degrees of freedom present 
in the thermal bath. 

Eq. (|273| ) summarizes all the ingredients we need to produce a baryon asymmetry in 
electroweak baryogenesis: 1) (the integral of ) a CP violating source A, 2) baryon number 
violation provided by the sphaleron transitions with rate r sp and 3) out-of-equilibrium 
conditions provided by the expanding bubble wall. 



7.2.3 Out-of-equilibrium field theory with a broad brush 

The next step in the computation of the baryon asymmetry is the evaluation of the CP 
violating sources for the right-handed stop number and the higgsino number. 

Non-equilibrium Quantum Field Theory provides us with the necessary tools to write 
down a set of quantum Boltzmann equations (QBE's) describing the local particle densities 
and automatically incorporating the CP violating sources. The most appropriate extension 
of the field theory to deal with these issues is to generalize the time contour of integration to a 
closed time-path (CTP). The CTP formalism is a powerful Green's function formulation for 
describing non-equilibrium phenomena in field theory, it leads to a complete non-equilibrium 
quantum kinetic theory approach and to a rigorous computation of the CP violating sources 
for the stop and the Higgsino numbers |109| , |110| , |lllf| . What is more relevant, though, is that 
the CP violating sources- and more generally the particle number changing interactions- 
built up from the CTP formalism are characterized by "memory" effects which are typical 
of the quantum transport theory |32] , [5S[ ]. CP violating sources are built up when right- 
handed stops and Higgsinos scatter off the advancing Higgs bubble wall and CP is violated 
at the vertices of interactions. In the classical kinetic theory the "scattering term" does 
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Figure 12: The contour C in the (Re t, Im t) plane proper to the CTP formalism 



not include any integral over the past history of the system. This is equivalent to assuming 
that any collision in the plasma does not depend upon the previous ones. On the contrary, 
the quantum approach reveals that the CP violating source is manifestly non-Mar kovian. 

We will now briefly present some of the basic features of the non-equilibrium quantum 
field theory based on the Schwinger-Keldysh formulation [115, |6l|]. The interested reader is 
referred to the excellent review by Chou et si. [^] for a more comprehensive discussion. 

Since we need the temporal evolution of the particle asymmetries with definite initial 
conditions and not simply the transition amplitude of particle reactions, the ordinary equi- 
librium quantum field theory at finite temperature is not the appropriate tool. The most 
appropriate extension of the field theory to deal with nonequilibrium phenomena amounts 
to generalize the time contour of integration to a closed-time path. More precisely, the time 
integration contour is deformed to run from — oo to +00 and back to —00, see Fig. 12. 

The CTP formalism (often dubbed as in-in formalism) is a powerful Green's function 
formulation for describing non-equilibrium phenomena in field theory. It allows to describe 
phase-transition phenomena and to obtain a self-consistent set of quantum Boltzmann equa- 
tions. The formalism yields various quantum averages of operators evaluated in the in-state 
without specifying the out-state. On the contrary, the ordinary quantum field theory (often 
dubbed as in-out formalism) yields quantum averages of the operators evaluated with an 
in-state at one end and an out-state at the other. 

Because of the time contour deformation, the partition function in the in-in formalism 
for a complex scalar field is defined to be 



Z 



J, J f 



Tr 



T 



exp 
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Tr 



x T_ exp 



T+ exp 



+ J 



t At 



P 



(274) 



where the sufhc C in the integral denotes that the time integration contour runs from minus 
infinity to plus infinity and then back to minus infinity again. The symbol p represents the 
initial density matrix and the fields are in the Heisenberg picture and defined on this closed 
time contour. As with the Euclidean time formulation, scalar (fermionic) fields <p are still 
periodic (anti-periodic) in time, but with <f>(t, x) = 4>{t — if3, x), (3 = 1/T. The temperature 
appears due to boundary condition, but time is now explicitly present in the integration 
contour. 

We must now identify field variables with arguments on the positive or negative di- 
rectional branches of the time path. This doubling of field variables leads to six different 
real-time propagators on the contour p2|| . These six propagators are not independent, but 
using all of them simplifies the notation. For a generic bosonic charged scalar field cj> they 
are defined as 



Gl (x,y) 
(x,y) 

G\{x,y) 
G r 4x,y) 



-i(<f>(x)<f>\y)}, 
-i(<j>\y)<j>{x)), 

9(x,y)G^(x,y) + 9(y,x)G^(x,y), 
9(y,x)G^(x,y) + 9(x,y)G$(x,y), 
G% — G^ = G^ — G^, G^(x, y) = G^ — G^ 



G< - Gt 



(275) 



where the last two Green functions are the retarded and advanced Green functions respec- 
tively and 9(x,y) = 9(t x — t y ) is the step function. For a generic fermion field ip the six 
different propagators are analogously defined as 



G% (x,y) 
G> (x,y) 
G%{x,y) 
G%{x,y) 
G r ^(x,y) 



-i(ip(x)ip(y)), 
+i$(y)4>{x)), 

9(x,y)G^(x,y)+9(y,x)G^(x,y), 
9(y,x)G^(x,y)+9(x,y)G^(x,y), 

G%{x,y)-- 



G^ 



G^, 



G^p, 



G 



G^p 



G 
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For equilibrium phenomena, the brackets (• • •) imply a thermodynamic average over all the 
possible states of the system. While for homogeneous systems in equilibrium, the Green 
functions depend only upon the difference of their arguments (x, y) = (x — y) and there 
is no dependence upon (x + y), for systems out of equilibrium, the definitions ((275) and 
(|276| ) have a different meaning. The concept of thermodynamic averaging is now ill-defined. 
Instead, the bracket means the need to average over all the available states of the system 
for the non-equilibrium distributions. Furthermore, the arguments of the Green functions 
(x,y) are not usually given as the difference (x — y). For example, non-equilibrium could 
be caused by transients which make the Green functions depend upon (t x ,t y ) rather than 

(t x — ty). 
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For interacting systems whether in equilibrium or not, one must define and calculate 
self-energy functions. Again, there are six of them: E , £*, S < , S > , E r and E a . The same 



relationships exist among them as for the Green functions in ( p75| ) and (276), such as 

S r = E* - E< = E> - E*, E a = E* - E> = E< - E*. (277) 

The self-energies are incorporated into the Green functions through the use of Dyson's 
equations. A useful notation may be introduced which expresses four of the six Green 



functions as the elements of two-by-two matrices [31] 



6 =l G> % ' S = \t % )' < 278 > 



where the upper signs refer to bosonic case and the lower signs to fermionic case. For 
systems either in equilibrium or non-equilibrium, Dyson's equation is most easily expressed 
by using the matrix notation 

G(x, y) = G°(x, y) + J d 4 x 3 J d 4 x 4 G°(x, x 3 )Z(x 3 ,x A )G(x A , y), (279) 

where the superscript "0" on the Green functions means to use those for noninteracting 
system. This equation appears quite formidable; however, some simple expressions may be 
obtained for the respective Green functions. It is useful to notice that Dyson's equation can 
be written in an alternate form, instead of ( |279| ), with G° on the right in the interaction 
terms, 

G(x, y) = G°(x, y) + J d*x 3 J d*x A G(x, x 3 )E(x 3 , x A )G°(x A , y). (280) 



Equations. ( |279|) and ( |280| ) are the starting points to derive the quantum Boltzmann equa- 
tions describing the temporal evolution of the CP violating particle density asymmetries. 

7.2.4 The quantum Boltzmann equations 

Our goal now is to find the QBE for the generic bosonic CP violating current 

(Jf(x)) = i(^(x) T x <Kx)) = [n^x), J^x)] . (281) 

The zero-component of this current represents the number density of particles minus 
the number density of antiparticles and is therefore the quantity which enters the diffusion 
equations of supersymmetric electroweak baryogenesis. 

Since the CP violating current can be expressed in terms of the Green function G^(x, y) 

as 

W*)> = - (5 - ^) G t^y)\ x=y > ( 282 ) 

the problem is reduced to find the QBE for the interacting Green function G^(x,y) when 
the system is not in equilibrium. This equation can be found from ( |279| ) by operating by 
(n x +m 2 ^J on both sides of the equation. Here m represents the bare mass term of the field 

(j). On the right-hand side, this operator acts only on G^ 

(a x +m 2 ) G^x, y) = 5^ (x, y)h + / d^^x, x 3 )G^(x 3 , y), (283) 



71 



where / is the identity matrix. It is useful to also have an equation of motion for the other 
variable y. This is obtained from ( |280| ) by operating by (o y +m 2 ^ on both sides of the 
equation. We obtain 



G<p(x,y) (a y +m 2 ) = <5 (4) (x,y)/ 4 + / d 4 x 3 G^,(x, xs)T,a(x 3 , y) 



(284) 



The two equations fl283|) and ( |284| ) are the starting point for the derivation of the QBE for 
the particle asymmetries. Let us extract from ( |283| ) and (284) the equations of motions for 
the Green function G5(x,y) 



(n x +m 2 ) G<(x,y) 



d 4 x 3 [^(x,x 3 )G<(x 3 ,y) - Z<(x,x 3 )Gl(x 3 ,y)] , (285) 
d 4 x 3 [G^(x,x 3 )^<(x 3 , y ) - G<(x,x 3 )^l(x 3 ,y)] . (286) 

If we now substract the two equations and make the identification x = y, the left-hand side 
is given by 



(d%-d%) G<(x,y) 



x=y 



d_4iX) 



dn 4 
dT 



and the QBE for the particle density asymmetry is therefore obtained to be 
dn^X) 



dT 



-+ v -MX) 



d x 3 



A^A ~ 2->A, ^A ~ ^A^A ~ A ^A 



x=y 



where we have defined the centre-of-mass coordinate system 

X = (T,X) = hx + y), (t,r)=x-y. 



(287) 



(288) 



(289) 



Notice that T now means the centre-of-mass time and not temperature. The identification 
x = y in Eq. (|288j ) is therefore equivalent to require t = r = 0. 

In order to examine the "scattering term" on the right-hand side of Eq. ( |288| ) , the 
first step is to restore all the variable arguments. Setting x = y in the original notation of 
T, ( f > (x,x 3 )G ( f > (x 3 ,y) gives (X, x 3 )(x 3l X) for the pair of arguments 



dn^X) 
dT 



+ V -U{X) 



d A x 3 \YL (X,x 3 )G<(x 3 ,X) - Z<(X,x 3 )Gi(x 3 ,X) 



+ Gi(X, x 3 )S< (x 3 ,X) - G< (X, x 3 )4(* 3 , X) 



(290) 



The next step is to employ the definitions in ( [275D to express the time-ordered functions 
G^,, G^,, and in terms of G^, G^, and G^. Then the time integrals are separated 
into whether t 3 > T or t 3 < T and the right-hand side of Eq. (|290|) reads 



d 4 x 3 {0(T-t 3 ) 



V>/^< i /"<<V> V</^> /^<>V< 
^A ^A ~T ^A ^A ~ ^A ~~ ^A ^A 



+ 0{t 3 -T) 



S</^r< i /"i<v< V</^< /^*<V< 
A ^A l ^A ^A ~ ^A ^A ~ ^A ^A 



(291) 
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The term with t 3 > T all cancel, leaving T > t%. Rearranging these terms gives fl09| , [11 
dn^X) 



w 



T 

dT 

G>(X,x 3 )X<(x 3 ,X) + G<(X,x 3 )X>(x 3 ,X) -X<(X,x 3 )G>(x 3 ,X) . (292) 



-+V / d 3 x 3 f dt 3 \z>(X,x 3 )G<(x 3 ,X) 

J J —00 L 



This equation is the QBE for the particle density asymmetry and it can be explicitly checked 
that, in the particular case in which interactions conserve the number of particles and the 
latter are neither created nor destroyed, the number asymmetry is conserved and obeys 

the equation of continuity dn ( p/dT+ V -jf = 0. During the production of the baryon 
asymmetry, however, particle asymmetries are not conserved. This occurs because the 
interactions themselves do not conserve the particle number asymmetries and there is some 
source of CP violation in the system. The right-hand side of Eq. (|292j ) , through the general 
form of the self-energy Ea, contains all the information necessary to describe the temporal 
evolution of the particle density asymmetries: particle number changing reactions and CP 
violating source terms, which will pop out from the corresponding self-energy T,qp- If the 
interactions of the system do not violate CP, there will be no CP violating sources and the 
final baryon asymmetry produced during supersymmetric baryogenesis will be vanishing. 



The kinetic Eq. ( 292 ) has an obvious interpretation in terms of gain and loss processes. 
What is unusual, however, is the presence of the integral over the time: the equation 
is manifestly non-Markovian. Only the assumption that the relaxation time scale of the 
particle asymmetry is much longer than the time scale of the non-local kernels leads to a 
Markovian description. A further approximation, i.e. taking the upper limit of the time 
integral to T — ► 00, leads to the familiar Boltzmann equation. The physical interpretation 
of the integral over the past history of the system is straightforward: it leads to the typical 
"memory" effects which are observed in quantum transport theory |58| . In the classical 
kinetic theory the "scattering term" does not include any integral over the past history of 
the system which is equivalent to assume that any collision in the plasma does not depend 
upon the previous ones. On the contrary, quantum distributions posses strong memory 
effects and the thermalization rate obtained from the quantum transport theory may be 
substantially longer than the one obtained from the classical kinetic theory. As shown in 



[ IPS , |l 10| , |lll| , memory effects play a fundamental role in the determination of the CP 



violating sources which fuel baryogenesis when transport properties allow the CP violating 
charges to diffuse in front of the bubble wall separating the broken from the unbroken phase 
at the electroweak phase transition. 

Notice that so far we have not made any approximation and the computation is therefore 
valid for all shapes and sizes of the bubble wall expanding in the thermal bath during a 
first-order electroweak phase transition. 

Let us now focus on the generic fermionic CP violating current. It reads 

(Jf(x)) = $(x)yV(z)> = [n^x), Xp(x)] , (293) 

where i/j indicates a Dirac fermion and 7^ represent the usual Dirac matrices. Again, 
the zero-component of this current represents the number density of particles minus 
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the number density of antiparticles and is therefore the relevant quantity for the diffusion 
equations of supersymmetric electroweak baryogenesis. 

We want to find a couple of equations of motion for the interacting fermionic Green 
function G 1 p(x,y) when the system is not in equilibrium. Such equations may be found 

by applying the operators I i jd x —M I and I i /) y +M I on both sides of Eqs. ( p79| ) and 



( |28C| ) , respectively. Here M represents the bare mass term of the fermion tp. We find 

f i%-M]G^{x,y) = 6^(x,y)h+ [ d 4 x 3 ^(x, x 3 )G^(x 3 ,y), (294) 

(295) 



5 (4) (x,y)7 4 - / d 4 x 3 G^(x,x 3 )T,^(x 3 ,y). 



G^(x,y) U jd y +M 

We can now take the trace over the spinorial indeces of both sides of the equations, sum up 
the two equations above and finally extract the equation of motion for the Green function 



Tr 



d A x 3 Tr Y^(x,x 3 )G\p(x 3 ,y) - T^(x,x 3 )G^{x 3 ,y) 



G ^( x , x 3)^(x 3 ,y) + G \i ) {x,x 3 )T^(x 3 ,y) 



(296) 



Making use of the centre-of-mass coordinate system, we can work out the left-hand side of 
Eq. (|296|) 



Tr 



l x G>(T,X,t, r) + G>(T,X,t, r)i jd y 



t=r=0 



-^$(*)YV(*)> 
d 

dx» u < p 



t=r=0 



Jtt. 



(297) 



The next step is to employ the definitions in fl276| ) to express the time-ordered functions 
G^p, G^, E^,, and in terms of G5, G^, £^ and G^. The computation goes along the 
same lines as the analysis made in the previous section and we get [109, 110, 111| 



dT 



v -M x ) 



d 3 x 3 



dU Tr 



T>{X,x 3 )G<{x 3 ,X) 



G>(X,x 3 )Z<(x 3 ,X) + G<(X,x 3 )Z>(x 3 ,X)-i:<(X,x 3 )G>(x 3 ,X)\. (298) 

This is the "diffusion" equation describing the temporal evolution of a generic fermionic 
number asymmetry rw,. As for the bosonic case, all the information regarding particle 
number violating interactions and CP violating sources are stored in the self-energy E^,. 



7.2.5 The CP violating source for higgsinos and the final baryon asymmetry 

As we mentioned, a strongly first order electroweak phase transition can be achieved in 
the presence of a top squark lighter than the top quark. In order to naturally suppress its 
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contribution to the parameter Ap and hence preserve a good agreement with the precision 
measurements at LEP, it should be mainly right-handed. This can be achieved if the left- 
handed stop soft supersymmetry breaking mass m,Q is much larger than Mz- Under this 
assumption, however, the right-handed stop contribution to the baryon asymmetry results 
to be negligible. We will concentrate, therefore, only on the CP violating source for the 
Higgsino. 

The Higgs fermion current associated with neutral and charged Higgsinos can be written 

as 

Jt = H'fH (299) 




H = (300) 



where H is the Dirac spinor 



and H 2 = H® {H}), Hi = H± (-Hf ) for neutral (charged) Higgsinos. The processes in 
the plasma which change the Higgsino number are the ones induced by the top Yukawa 
coupling and by interactions with the Higgs profile. The interactions among the charginos 
and the charged Higgsinos which are responsible for the CP violating source in the diffusion 
equation for the Higgs fermion number read 

C = -g 2 {H [ Vl {x)P L + e id »v 2 (x)P R ] w) + h.c, (301) 

where 9^ is the phase of the /^-parameter and we have indicated {Hf(x)} by Vi(x), i = 1,2. 
Analogously, the interactions among the Bino, the H^-ino and the neutral Higgsinos are 

C = -\ {~H° [vi{x)P L + e^v 2 (x)P R ] (g 2 W 3 - 9i b) } + h.c. (302) 

To compute the source for the Higgs fermion number 7- we perform a "Higgs insertion 
expansion" around the symmetric phase. At the lowest level of perturbation, the inter- 
actions of the charged Higgsino induce a contribution to the self-energy of the form (and 

CP > \ 

analogously for the other component oS~ ' ) 



«5E| P ' < (x, y) = gS P (x, y)P L G°~ < (x, y)P L + g* P (x, y)P R G^ < (x, y)P R , (303) 



where 



9cp( x iV) = 9 2 2 vi{x)v 2 {y)e id », 

gg P (x,y) = dv^yfaWe*". (304) 

We have approximated the exact Green function of winos by the equilibrium Green 
function in the unbroken phase G . This is because any departure from thermal equilibrium 
distribution functions is caused at a given point by the passage of the wall and, therefore, 
is 0{vj). Since we will show that the source is already linear in v u , working with thermal 
equilibrium Green functions in the unbroken phase amounts to ignoring terms of higher 
order in v u . This is accurate as long as the bubble wall is moving slowly in the plasma. 
Similar formulae hold for the neutral Higgsinos. 
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The dispersion relations of charginos and neutralinos are changed by high temperature 
corrections |123( |. Even though fermionic dispersion relations are highly nontrivial, especially 
when dealing with Majorana fermions |107| 1, relatively simple expressions for the equilibrium 
fermionic spectral functions may be given in the limit in which the damping rate is smaller 
than the typical self-energy of the fermionic excitation . If we now insert the expressions 
(BOg) and (gOp into the QBE (Egg), we get the CP violating source EM pM 



7 



d 6 x 3 



dU Tr 



<5S| p ' > (X,x 3 )Gj < ( 3 ;3,X) - G^{X,x 3 )5T^^{x 3 ,X) 



y0,>. 



^CP,<, 



lf(X,x 3 )6i:f>(x 3 ,X)-6Zf>< 
which contains in the integrand the following function 



+ G»f(X,x 3 )6XY'>(x 3 ,X) - 6Z^(X,x 3 )Gf(x 3 ,X) 



(305) 



g C p(X,x 3 )+ggp(X,x 3 )-g C p(x 3 ,X)-gg P (x 3 ,X) = 2is^e^ [v 2 (X)vi(x 3 ) - Vi(X)v 2 (x 3 )] , 

(306) 

which vanishes if Im(/i) = and if the tan/3(x) is a constant along the Higgs profile. 

In order to deal with analytic expressions, we can work out the thick wall limit and 
simplify the expressions obtained above by performing a derivative expansion 



H(x 3 ) = 



n=0 



1 



n\ d{X^)'' 



MX) K-iT. 



(307) 



The term with no derivatives vanishes in the expansion ( 307 ), v 2 (X)vi(X) — vi(X)v2(X) = 
0, which means that the static term in the derivative expansion does not contribute to the 
source. For a smooth Higgs profile, the derivatives with respect to the time coordinate and 
n > 1 are associated with higher powers of v^/L^, where v w and are the velocity and 
the width of the bubble wall, respectively. Since the typical time scale of the processes 
giving rise to the source is given by the thermalization time of the higgsinos 1/T~, the 
approximation is good for values of L^T-fv^ 3> 1. In other words, this expansion is valid 
only when the mean free path of the higgsinos in the plasma is smaller than the scale 
of variation of the Higgs background determined by the wall thickness, L u , and the wall 



velocity v u . The term corresponding to n = 1 in the expansion ( 307 ) gives a contribution 
to the source proportional to the function 



v 1 (X)d^v 2 (X) - v 2 (X)d» Vl (X) = v\X)d» x P(X) 



(308) 



which should vanish smoothly for values of X outside the bubble wall. Here we have denoted 
v 2 = v\ + v\. Since the variation of the Higgs fields is due to the expansion of the bubble 
wall through the thermal bath, the source 7^ will be linear in v u . The corresponding 
contribution tot he CP violating source reads 



75 (X)=Im(/i) v\X)f}(X) 



3M 2 g\ Jf + Mi g\ Zf 



(309) 



where 



X 



!'OG 



dk- 



k 2 



r .... 

Jo 2vr 2 o;~a;~ 
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(l - 2Re(f~)) J( W ~,r~,u,~,r~) + (l - 2Re(/|)) , F~ , w ~, T~) 
+ 2(lm(/|)+Im(4))G(^,r~, W ~,r~) 



and w 2 - ~ = k 2 + |u| 2 (Mf) while /2. ~ = 1/ [exp (u~,~JT + iT~,~. IT) 

H(W) lrl V 2 ' J H(W) ' L tf(Wy tf(Wy / 

functions / and G are given by 



(310) 



The 



I(a,b,c,d) = 



2[(a + C y + (b + d) 2 } 



1 



G(a, 6, c, c/) 



2[(a-c) 2 + (6 + d) 2 ; 
1 1 
2[(a + c) 2 + (6 + d) 2 ; 
1 1 
2[(a-c) 2 + (6 + d) 2 ; 



sm 



sm 



cos 



cos 



2arctan 
2arctan 
2arctan 
2arctan 



a + c 



b + d 
a — c 



b + d_ 
a + c 
6 + d 
a — c 
6 + d 



(311) 



uj~. This resonant be- 
W 



Notice that the function G{oj~, F~, uj^, r^p.) has a peak for w~ 

haviour is associated to the fact that the Higgs background is carrying a very low momentum 
(of order of the inverse of the bubble wall width L w ) and to the possibility of absorption 
or emission of Higgs quanta by the propagating supersymmetric particles. The resonance 
can only take place when the higgsino and the wino do not differ too much in mass. By 
using the Uncertainty Principle, it is easy to understand that the width of this resonance 
is expected to be proportional to the thermalization rate of the particles giving rise to the 
baryon asymmetry. 

The damping rate of charged and neutral Higgsinos is expected to be of the order of 
5 x 1CT 2 T. The Bino contribution may be obtained from the above expressions by replacing 
M2 by M\. The CP violating source for the Higgs fermion number is enhanced if M2, M\ ~ 
[i and low momentum particles are transmitted over the distance L^. This means that 
the classical approximation is not entirely adequate to describe the quantum interference 
nature of CP violation and only a quantum approach is suitable for the computation of the 
building up of the CP violating sources. Notice that the source is built up integrating over 
all the history of the system. This leads to "memory effect" that are responsible for some 
enhancement of the final baryon asymmetry. These memory effects lead to "relaxation" 
times for the CP violating sources which are typically longer than the ones dictated by 
the thermalization rates of the particles in the thermal bath. In fact, this observation is 
valid for all the processes described by the "scattering" term in the right-handed side of 
the quantum diffusion equations. The slowdown of the relaxation processes may help to 
keep the system out of equilibrium for longer times and therefore enhance the final baryon 
asymmetry. There are two more reasons why one should expect quantum relaxation times to 
be longer than the ones predicted by the classical approach. First, the decay of the Green's 
functions as functions of the difference of the time arguments: an exponential decay is 
found in thermal equilibrium when one ignore the frequency dependence of self-energies in 
the spectral functions, e.g. \G > (k,t,t')\ ~ |G > (k)| x exp [— T(k, uj)\t — t'\]. The decay of the 
Green's functions restrict the range of the time integration for the scattering term, reduces 
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the integrals and, therefore, the change of the local particle number densities as a function 
of time. The second effect is the rather different oscillatory behaviour of the functions G > 
and G < for a given momentum, as functions of the time argument difference. 

As we have previously mentioned, the final baryon asymmetry fl273| ) depends sensitively 
on the parameter A. The parameter A computed from the higgsino source is 

A oc J JL tJ H I, 
D A+ 







I = / (f!i!)1n)-^e- A +V duv 2 (u)—^, (312) 
du Jo du 



where /(fcj) is a coefficient depending upon the number of degrees of freedom present in 
the thermal bath. The integral I has been computed including two-loop effects in ref. [106] 



and results to be / ~ 10 2 for rriA = 150-200 GeV. The final baryon asymmetry turns out 
to be gig 

^^^M) 4xl0 _n (313) 

for v w ~ 1. It is intriguing that these small values of the phases are perfectly consistent 
with the constraints from the electric dipole moment of the neutron and squarks of the first 
and second generation as light as ~ 100 GeV may be tolerated. 



8 Conclusions 

In these lectures we have learned that cosmology provides really strong arguments in favour 
of the nonconservation of the baryon number. The SM of weak interactions, which is so 
successfull in explaining the experimental data obtained at accelerator machines operating 
at energy scales of about 100 GeV, seems unable to explain the observed baryon asymmetry 
of the Universe. This is a very strong indication that there is some new, yet undiscovered, 
physics beyond the SM. We do not know whether this is just the low energy supersymmetric 
extension of the SM. If so, we can draw tight constraints on the Higgs spectrum of the 
MSSM and the next generation of accelerator machines, such as LHC, will tell us if this is 
a tenable option. It might be that this cosmological puzzle has been taken care of by some 
new physics at energy scales much higher than the weak scale, the GUT scale, as suggested 
by gauge coupling unification. Even though this option is not testable at particle colliders, 
the most striking evidence of baryon number violation might come from the detection of 
proton decay. It is very exciting that in the next few years we will be able to confirm (or 
disprove) some of the theories of baryogenesis. 
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Solution to the Exercises 



1) Suppose that the baryon number B is conserved by the interactions. This means 
that the baryon number commutes with the hamiltonian of the system H, [B,H] = 0. 
Therefore, supposing that B(to)=0, we have B(t) oc / t * [B,H]dt' = at all times and no 
baryon number production may take place. 



2) Let us consider a model universe with two components: inflaton field energy, p^ 
and radiation energy density, pr. We will assume that the decay rate of the inflaton field 
energy density is T^. We will also assume that the light degrees of freedom are in local 
thermodynamic equilibrium. 

With the above assumptions, the Boltzmann equations describing the redshift and in- 
terchange in the energy density among the different components is 

p^ + zup^ + = o 

p R + AH p R - I> = 0, (314) 



where dot denotes time derivative. 

It is useful to introduce the dimensionless constant, a<j> defined in terms of as 

1^ = c^ilV (315) 

For a reheat temperature much smaller than M^, must be small. 

It is also convenient to work with dimensionless quantities that can absorb the effect of 
expansion of the universe. This may be accomplished with the definitions 

$ = p^M^a 3 ; R = p R a\ (316) 

It is also convenient to use the scale factor, rather than time, for the independent variable, 
so we define a variable x = aM^. With this choice the system of equations can be written 
as (prime denotes d/dx) 

$ = -Cl ; $ 
y/$X + R 

R' = ci ; X (317) 
V$x + R 



The constant c\ is given by 



cl = v^M: a ^ (318) 



It is straightforward to solve the system of equations in Eq. (317) with initial conditions 
at x = xi of R{xi) = and $(xi) = It is convenient to express p${x = xf) in terms of 
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the expansion rate at xi, which leads to 



^1 Hi , 



8vrM 2 M 2 ' 



(319) 



The numerical value of xi is irrelevant. 

Before solving the system of equations, it is useful to consider the early-time solution for 
R. Here, by early time, we mean H 3> i.e., before a significant fraction of the comoving 
coherent energy density is converted to radation. At early times 3> ~ and R ~ X ~ 0, 
so the equation for R' becomes R' = cix 3 ^ 2 ^ 1 / 2 . Thus, the early time solution for R is 
simple to obtain: 



R ~ - Cl (; 



5/2 _ X 5/2N, a l/2 



(# » 



Now we may express T in terms of i? to yield the early-time solution for T: 



T 




x\ J 



-3/2 



1/4 



(320) 



(321) 



Thus, T has a maximum value of 

Tmax 

Mrh 



0.77 



0.77a 



which is obtained at x/xi = (8/3) 2//5 
Trh and obtain 

Tmax 




2tt 3 9*J \ Ml 
1.48. It is also possible to express 



(322) 



in terms of 



0.77 



/ 9 \ 1 ^(H I Mj>\ 



1/4 



(323) 



Trh ""' \2ir 3 gJ \Trh ) 

For an illustration, in the simplest model of chaotic inflation Hj ~ MpMj, with ~ 
10 13 GeV, which leads to Tmax / Trh ~ 2 x 10 3 (200/5*) 1 / 4 for Tr H = 10 9 GeV. 

We can see from Eq. ( p20[ ) that for x/x/ > 1, in the early-time regime T scales as a~ 3 / 8 . 
So entropy is created in the early-time regime. So if one is producing a massive particle 
during reheating it is necessary to take into account the fact that the maximum temperature 
is greater than Trh, and during the early-time evolution T oc a~ 3 / 8 . 



3) The euclidean action is given by 

f R 

S3 = 4-7T / r 2 dr 
Jo 



dr 



+ v{4>) 



(324) 



where R is the radius of the bubble. If we now indicate by 5R the thickness of the bubble 
wall and AV = V{4> 2 ) - V{4>\) < 0, we get 



S 3 ~ 2irR 2 5R 



5RJ 



2 rn 4irR 3 AV 
5R + , 



(325) 
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where 8<f> = (f>2 — <pi- Suppose now that the bubbles are thick, that is 5R ~ R. In such a 
case 

5 3 ~ 2vr R(5ct>) 2 + 4?ri ^ AF . (326) 
The critical radius R c is obtained as the maximum of the action ( |326| ) 

R c ~ -^=. (327) 



4 ) We introduce a chemical potential for any particle which takes part to fast processes, 
and then reduce the number of linearly independent chemical potentials by solving the 
corresponding system of equations. Finally, we can express the abundances of any particle in 
equilibrium in terms of the remaining linear independent chemical potentials, corresponding 
to the conserved charges of the system. 

Since strong interactions are in equilibrium inside the bubble wall, we can chose the 
same chemical potential for quarks of the same flavour but different color, and set to zero 
the chemical potential for gluons. Moreover, since inside the bubble wall SU{2)i x U(1)y 
is broken, the chemical potential for the neutral Higgs scalars vanishes^. 

The other fast processes, and the corresponding chemical potential equations are: 

i) top Yukawa: 

t L + H$ ^t R + g, (Vt L = MtJ, / 328 n 

b L + H + ^t R + g, (fi tR = fi bL +fi H +), 



ii) SU(2)l flavour diagonal: 



(/V r — Mel + 



L 



'L 

< ~ 4 + W+, (m** = M 4 + A**"-), {i = ! 2 3) . (329) 

H°~H-+W+, {ijl h -= -im w+ ), 

Neutral current gauge interactions are also in equilibrium, so we have zero chemical potential 
for the photon and the Z boson. 

Imposing the above constraints, we can reduce the number of independent chemical po- 
tentials to four, fit L , Hu L = 1/2 5Zi=i AV ; an d He L = ^/3J2i=i AV • These quantities 
correspond to the four linearly independent conserved charges of the system. Choosing the 
basis Q, (B — L), (B + L), and BP = B3 — \j2{B\ +B2), where the primes indicate that only 
particles in equilibrium contribute to the various charges, and introducing the respective 
chemical potentials, we can go to the new basis using the relations 




3fi tL + 2fi UL - 3fi eL + 11 fj, w +, 

(330) 



3^i L + Afi UL - 6fi eL - 6fi w + , 



3nt L + 4:fi UL + 6/i eL , 
3^ L - 2fi UL . 



4 



This is true if chirality flip interactions, or processes like Z —> Z*h, are sufficiently fast. 
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If sphaleron transitions were fast, then we could eliminate a further chemical potential 
through the constraint 

3 3 3 

i=l i=l i=l 

In this case, the value of {B + L) would be determined by that of the other three charges 
according to the relation 

S 7 i q 
(B + L) E o = — Q H BP (B-L). (332) 

The above result should not come as a surprise, since we already know that a non zero value 
for B — L gives rise to a non zero (B + L) at equilibrium. Stated in other words, sphaleron 
transitions erase the baryon asymmetry only if any conserved charge of the system has 
vanishing thermal average, otherwise the equilibrium point lies at (B + L)eq +~ 0. 
At high temperature (/tj <C T) the free energy of the system is given by 



F -Tl 

r ~ 12 



3 /4 + 3 /4 + 6 /4 + M L + M R + 6/4 + 3/4 



+6/4+ + 2/4+ + 2/4 ? + 2/i^o 



Using ( p28| ), (|329[ ) and ( p30| ) to express the chemical potentials in terms of the four conserved 
charges in ( |330| ) we obtain the free energy as a function of the density of (B + L) = 
fi B+L T 2 /6, 

F [(B + L)\ = 0.46 [(jB + L) ~ T ( f + L)eq] + constant terms, (334) 

where the "constant terms" depend on Q, {B — L), and BP but not on (B + L), and 
(B + L) E q is given by ( |332| ). 

From the above expression, we can see that 



5) The contribution to the E parameter of the potential ( 231 ) is given generically by, 
see Eq. (Egg), 



— E 

127T 4" 



mf^+n^T) 



3/2 



(336) 



for a generic bosonic particle i with plasma mass Ilj(r). In the case of the right-handed 
stop, we get 



5E 



-2 N r . 



12vr 



m 2 ~+ U R (T) 



3/2 



(337) 



where A^ c = 3 is the number of color and m~ is given in ( p51| ). The upper bound on the 
contribution to the E parameter from the right-handed stops is obtained when mfj < and 



m 



off 



mfj + U R (T) ~ 0. This gives 



h 3 t 



3/2 



(338) 
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Using now the fact that mt 



h t v and that (<f>(T c ))/T c = 2E/X ~ 4:V 2 E/m 2 h , we get Eq. 
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